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SUMMARY 


This  research  report  describes  studies  done  by  C-E-I-R,  Inc., 
under  Nonr  3333 (00) ,  which  are  concerned  with  the  impact  of  obsolescence 
on  inventories  of  Naval  supplies.  Obsolescence  is  examined  with  the 
objective  in  mind  of  defining  it  operationally.  Several  types  of 
dynamic  models  are  developed,  making  use  of  the  operationally  defined 
motion.  Some  of  these  models  employ  Bayesian  procedures  for  assessment 
of  demand  rates.  A  statistical  analysis  of  some  Navy  supply  data, 
which  was  carried  out  in  an  attempt  to  cover  empirical  relationships 
between  obsolescence  and  other  factors,  is  described.  Finally  some 
changes  in  inventory  control  procedures  are  suggested. 


TABLE  OF  CONTENTS 


Chapter  I: 
Chapter  II: 

Chapter  III: 

Chapter  IV: 
Chapter  V: 

Chapter  VI: 
Chapter  VII: 

Chapter  VIII 
Appendices: 

References 


Introduction  1 

Obsolescence:  Definition  and 

Importance  3 

Survey  of  Inventory  Models  with  Explicit 
Treatment  of  Obsolescence  13 

A  Dynamic  Procurement-Disposal  Model  33 

Reproducibility  of  Probability 

Distributions  as  Related  to  Bayesian 

Inventory  Models  44 

Inventory  Models  with  Markov  Demand  56 

Some  Results  of  a  Statistical  Study  of 
Obsolescence  Data  77 

Summary  and  Conclusions  94 

A.  Policies  Resulting  from  Dynamic 

Procurement  -  Disposal  Model  Runs  A- 1 

B.  FORTRAN  Programs  for  Inventory 

Models  Developed  in  Chapters  4  and  6  B-l 

C.  Discussion  of  ASO  Procedures  C-l 


CHAPTER  I:  INTRODUCTION 


In  January  1961,  C-E-I-R,  Inc.,  was  authorized  by  the  Bureau  of 
Supplies  and  Accounts  to  undertake  a  study  of  the  economic  impact  of 
obsolescence  on  inventory  costs  and  control.  This  document  is  a  re¬ 
port  of  results. 

The  study  has  led  to  the  development  of  a  sequence  of  models  of 
inventory  costs,  computing  techniques  associated  with  these  models, 
and  a  proposal  for  installation  of  inventory  control  procedures  based 
on  the  models  described  herein.  The  models  deal  with  inventory  opera 
t ion  costs  under  a  variety  of  assumptions  as  to:  lead  time,  inven¬ 
tory  carrying  charges,  ordering  costs,  disposal  costs,  shortage 
costs,  probability  distribution  of  demand  for  an  item,  probability 
of  incidence  of  obsolescence  during  any  time  period  and  linkage  be¬ 
tween  demands  for  two  or  more  items. 

These  models  have  one  feature  in  common:  all  are  dynamic  models 
That  is,  a  steady  state  is  not  assumed;  the  possibility  exists  for 
variation  of  demand  and  obsolescence  distributions,  as  well  as  all 
the  costs  of  operation  of  the  inventory.  Several  of  the  more  sophis¬ 
ticated  models  incorporate  Bayesian  features.  Taken  together,  these 
two  features,  the  dynamic  program  and  the  Bayesian  features,  lead 
to  an  ability  to  adapt  in  a  continuous  fashion  to  changes  both  in 
the  structure  of  the  environment  and  in  a  priori  statements  whose 
validity  is  revealed  only  slowly  as  experience  is  generated. 
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Central  to  this  study  is  an  attempt  to  understand  and  model  the 
nature  and  impact  of  obsolescence  of  items  in  inventory.  The  follow¬ 
ing  chapter  is  concerned  with  a  discussion  of  what  is  meant  by  obso¬ 
lescence,  as  well  as  some  concern  with  the  magnitude  of  obsolescence 
of  inventories  of  Naval  Supplies. 
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CHAPTER  II:  OBSOLESCENCE:  DEFINITION "AND  'IMPORTANCE 


It  has  been  standard  practice  in  the  U.S.  Navy  supply  system  to 
deal  with  obsolescence  as  a  blanket  charge  to  be  assessed  equally 
across  all  Navy  Supply  items.  The  rate  at  which  this  charge  has 
customarily  been  assessed  is  10%  of  the  total  value  of  the  inven¬ 
tory  per  year. 

The  total  value  of  inventories  of  Naval  Aircraft  spare  parts 
and  supplies  held  by  the  Navy  amounts  to  approximately  $2.5  billion. 
Clearly,  the  total  value  of  all  inventories  held  by  the  Navy  is 
muck  larger  than  this  amount.  Consequently,  the  total  obsolescence 
charges  under  the  present  system  amount  to  several  hjindred  million 
dollars  per  year. 

It  seems  indisputable  that  better  understanding  of  the  specific 
impact  of  obsolescence,  and  the  development  of  inventory  control 
procedures  which  take  the  obsolescence  process  into  account  in  a 
realistic  fashion,  should  lead  to  considerable  economies  in  the  use 
of  Navy  Supply  budgets. 

Defining  the  term  obsolescence  is  not  quite  as  simple,  in  the 
context  of  a  study  of  the  sort  under  discussion,  as  it  would  be  for 
the  purposes  of  ordinary  discourse.  Something  is  obsolescent ,  ac¬ 
cording  to  Webster's  New  Collegiate  Dictionary  (1960  Edition),  when 
it  is  going  out  of  use.  In  ordinary  discourse,  it  is  assumed  that 
there  is  no  serious  problem  involved  in  ascertaining  whether  or  not 


something  is  going  out  of  use.  However,  for  our  purposes,  it  turns 
out  that  this  is  a  most  difficult  matter  to  settle  until  after  the 
fact . 

Establishing  that  something  is  obsolete,  that  is,  establishing 
that  something  has  alreadj  gone  out  of  use;  is  a  much  simpler  problem 

In  this  case,  there  are  several  circumstances  which  might  have 
led  to  such  a  situation. 

1.  An  item  has  gone  out  of  use  because  the  function  served 
by  that  item  is  no  longer  required;  i.e.;  demand  has  dis¬ 
appeared. 

2.  An  item  is  going  out  of  demand  because  whenever  a  unit  of 
the  item  wears  out,  breaks,  or  is  consumed,  the  item  is 
replaced  by  a  different  item  which  performs  similar  or 
identical  functions;  i.e.,  the  item  is  going  out  of  use. 

3.  An  item  has  gone  out  of  demand  because,  regardless  of  con¬ 
sumption,  wearout  or  breakage,  as  soon  as  replacements  are 
available,  all  units  of  that  item  currently  in  use  are  re¬ 
placed  by  units  of  other  items  which  serve  similar  or  iden¬ 
tical  functions,  that  is,  the  item  has  gone  out  of  use. 

So  far  as  the  supply  system  is  concerned,  being  out  of  use,  and  not 
being  in  demand  are  the  same  condition. 

Once  something  has  become  obsolete,  the  fact  is  clear.  But  ther 
is  not  much  which  can  be  done  at  such  a  time  to  face  the  fact,  beyond 
clearing  out  and  disposing  of  any  remaining  inventories  of  such  items 


something  is  going  out  of  use.  However,  for  our  purposes,  it  turns 
out  that  this  is  a  most  difficult  matter  to  settle  until  affer  the 
fact . 

Establishing  that  something  is  obsolete,  that  is,  establishing 
that  something  has  already  gone  out  of  use;  is  a  much  simpler  problem. 

In  this  case,  there  are  several  circumstances  which  might  have 
led  to  such  a  situation. 

1.  An  item  has  gone  out  of  use  because  the  function  served 
by  that  item  is  no  longer  required;  i.e.,  demand  has  dis¬ 
appeared. 

2.  An  item  is  going  out  of  demand  because  whenever  a  unit  of 
the  item  wears  out,  breaks,  or  is  consumed,  the  item  is 
replaced  by  a  different  item  which  performs  similar  or 
identical  functions;  i.e.,  the  item  is  going  out  of  use. 

3.  An  item  has  gone  out  of  demand  because,  regardless  of  con¬ 
sumption,  wearout  or  breakage,  as  soon  as  replacements  are 
available,  all  units  of  that  item  currently  in  use  are  re¬ 
placed  by  units  of  other  items  which  serve  similar  or  iden¬ 
tical  functions,  that  is,  the  item  has  gone  out  of  use. 

So  far  as  the  supply  system  is  concerned,  being  out  of  use,  and  not 
being  in  demand  are  the  same  condition. 

Once  something  has  become  obsolete,  the  fact  is  clear.  But  there 
is  not  much  which  can  be  done  at  such  a  time  to  face  the  fact,  beyond 
clearing  out  and  disposing  of  any  remaining  inventories  of  such  items. 
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However,  if  it  were  possible  to  detect  accurately  the  process  of  obso' 
lescence  in  its  earlier  stages,  while  it  has  not  run  its  course,  then 
there  might  be  some  steps  which  might  profitably  be  taken  by  an  in¬ 
ventory  manager.  Consider,  for  example,  the  following  very  simple 
case: 

Let: 

Carrying  Cost  =  h 
Stockout  Cost  =  h 
Liquidation  Cost  =  k 
Demand  =  D 

Probability  distribution  of  demand  =  p(D) 

P(D)  =  P(r)  =  Prob  (IKr) 

A.  It  is  not  known  that  obsolescence  will  occur  at  the  end  of  this 
period.  An  optimal  policy  is  followed,  then  obsolescence  occurs  and 
stock  is  liquidated. 

The  expected  inventory  costs  the  manager  can  plan  on  are: 

(2.1)  C  (I  )  =  h  Z°  (I  -D)  p (D)  +  it  L  (D-I  )  p  (D) 

0  D=0  0  D=I0+ 1  0 

If  I  is  optimal 
o 

(2.2)  C(Io  +  1)  -  C(Iq)  >  0 
C(IQ  -  1)  -  C(IQ)  >  0 


HeI 


n« 


(2.3) 
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(2.4)  C(I  +1)  -  h  1°  (I  -D>  p(D)  +  hj°  p(D)  +  R  S  ^  (D-I  )  p(D) 
o  D*0  o  D«0  D"I0+1  0 


'>?0+i p<B) 


(2.5)  0(Io+l)  -  0(Io)  -  P(I0)  (hMi)  -  «  >  0 


Similar ly, 


(2.6)  -  [c(Io-l)  -  C (Io) ]  -  P  (I0-l)  (hfn)  -  n  < 

<2’7> 


B.  It  is  known  that  obsolescence  will  occur  at  the  end  of  the 
period.  The  manager  plans  accordingly.  His  expected  costs  are: 

(2.8)  E(I*)  -hi*  (I*-D)  p(D)  +  x  (D-I*)  p(D)  +  kpl*  (I*-D)  p(D) 


(2.9)  E(I.)  -  (hfk)  Z  (I.-D)  p(D)  +n  _  Z  _  (D-I.)  p(D) 


D-0  * 


D-I#fl 


By  the  same  reasoning  as  (2.1)  to  (2.7) 


(2.io)  p(Vd  <  — ; <  p(i*) 


80  h  <  lo 


The  total  expected  costs  incurred  in  situation  A.  must  include 
liquidation  of  remaining  stock: 


(2.11)  E(I0)  ■  C(I0)  +  k  Z  «o-0)  p (D)  •  (hMO  L  (I0-D)  p(D) 


+  "  >&1  <d’i*)  p(d> 


(2.12)  E(Iq)  -  E(I*)  =  (hfk)  [jjEq  (I0-D)  p(D)-^0  (I*-D)  P  (D)J 

+*  [u=ijfi  *  D-lf+l<D'V  pO»] 

It  is  easily  shown  (see  Footnote  1)  that  (2.12)  is  positive. 


'’’Starting  with  (2.12):  adding  (hfk)  ^"^(I^D)  p(D)-J>=IZ+1(I^-D)p(D)| 

_  O  ic  J 


.-|+1  <V»  P0»] 


-  <W-fc)  D-A-l  <V“>  P(D)-  D-A-l  <VI»  P<») 

I  O  * 


we  get 

(2.13)  E(Io)-E(I*)=(hfk)jj0  (Io-D)  p(D)  -  J0  (VD)  p(D)j 

+  (^4it)|?=iL+i(D-Io)  P(D)-d=i2+1(D-I*)  pW] 

(2.14)  E(Io)-E(I*)=(hfk)(Io-I^)+(hfk+7t)  ^j£+1  (D-Iq)  P (D) 

-D=l”+1  <D-X*>  P<°>  '  D=Jfl  P(D)1 

o  *  J 

(2.15)  E(Io)-E(I^)=(hfk)(Io-I#)+(hfk-Ht)^D=I“+1  (I*-I0)  P(D) 

-d-4+1  <d’i*)  p<D)] 

(2.16)  «(Io)-E(I4)-(hPk)<Io-ip+0»*hr)  {(I*-1,)  [1_P,'Io)] 


A  (o-V  po)) 


-8- 


(2.17)  E(Io)-E(I*)  -  (hfk)(IQ-I#)  +  (hfk+n) 

[p<1*>  +  D-i^iP®]  -  o-jfli  OVpW) 

(2.18)  E(Io)-E(I*)  *  -  *(IQ-I*)  +  (hflrt*)(Io-I*)  P  (I*) 

T-o 

+  D-iJ+l  (V1*  ‘Dfl*>  P(D) 

(2.19)  E(Io)-E(I*)  =  -*  (Io-I*)  +  (hfkHt)(lo-l#)(  ) 

+  (hfk4*)  (Id-D)  p(D) 

(2.20)  E(I0)-E(I*)  =  -n  (IQ-I*)  +  (>*(I0-I*))  +  (hfk+n)  J°(I  -D)p(D) 

'— — —  _ /  V  ■L*~ 

>0 

We  have  made  use  of  the  fact  that 

P(I.)  >  *  and 

v  *  —  n+hhk 

Th:  If  I  >  J  and 

then  E(I)  >  E(J) 

We  now  consider  a  numerical  example: 
h  =  60 
n  =  500 
k  =  200 

Demand  is  Poisson  distributed,  with  X= 


>0 


I*<Io 


P(J)  > 


n+hfk 


-6  D 

6>P(D)  =  --  Dj6^ 
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«  500 

jt+h  "  560 


.8928  I0 


9 


it  _  500 
jr+h+k  860 


.5814  I*  «  6 


from  (2.18) 

9 

(2.18')  Ed^-EU*)  -  -500(3)  +  (860)  (3)  (.6063)  +  (806)  ^ ?  (9-D)  p  (D) 

=  -1500  +  (2580) (.6063)  +806  2  x  .1377  +  1  x  .1033 

=  369.4862 

Total  costs  under  situation  B  are  732.5640 

Thus,  there  is  clear  gain  of  the  order  of  507.  from  having  had  prior 
knowledge  of  the  impact  of  obsolescence. 

In  cases  where  items  go  out  of  demand  or  out  of  use  as  a  conse¬ 
quence  of  administrative  decisions,  it  might  be  easier  to  anticipate 
the  process  of  obsolescence  than  in  other  cases.  This  is  sc  because 
administrators  must  have  some  reasons  for  making  such  decisions,  and 
these  reasons  are  probably  open  to  examination.  However,  in  cases 
where  administrative  decisions  are  not  involved,  it  is  necessary  for 
inferences  to  be  made  about  the  future  state  of  demand  from  the  pres¬ 
ent  state  of  demand. 

Such  inferences  are  of  the  form:  given  that  demand  in  the  pres¬ 
ent  period  is  equal  to  y,  the  probability  that  demand  M  months 
from  now  will  be  less  than,  or  equal  to  y  is  P(y;y,M).  When  P  is 
sufficiently  high,  it  may  turn  out  that  the  cost  of  maintaining  the 
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item  in  inventory  for  demand  levels  as  low  as,  or  lower  than,  y  is 
more  than  the  expected  costs  which  would  be  incurred  if  the  system 
were  out  of  stock  of  the  item. 

In  such  a  circumstance  a  declaration  that  the  item  is  to  be 
treated  as  if  it  were  obsolescent  might  be  in  order.  This  would 
mean  taking  some  action  which  would  lead  to  one  of  the  three  conditions 
1) ,  2) ,  or  3) ,  which  are  described  above,  (on  page  4) . 

In  the  models  which  are  developed  and  discussed  in  the  follow¬ 
ing  sections  of  this  report  there  are,  basically,  two  approaches 
taken  to  the  notion  of  obsolescence. 

The  first,  a  simpler  treatment,  is  concerned  with  the  develop¬ 
ment  of  dynamic  disposal-procurement  policies  under  the  assumption  of  a 
probability  distribution  of  time  of  incidence  of  obsolescence  which  is 
known  a  priori,  and  which  is  not  changed  systematically  within  the 
mode 1 . 

This  model  is  very  flexible  so  far  as  input  data  are  concerned; 
any  form  of  obsolescence  probability  distribution,  demand  distribu¬ 
tion,  or  inventory  costs  may  be  used.  Consequently,  as  ancillary  in¬ 
formation  about  any  of  these  is  developed,  new  assumptions  can  be  in¬ 
troduced  into  the  model. 

Numerical  analysis  of  this  model  has  been  carried  out  fairly 
extensively.  The  results  of  the  numerical  analysis  tend  to  support 
conjectures  as  to  the  extensibility  and  greater  generality  of  some 
theorems  which  have  appeared  in  the  literature  recently.  Moreover, 


this  model  appears  promising  as  the  basis  for  a  first  step  toward  the 
rationalization  of  procurement-disposal  stock  control  of  Naval 
inventories. 

The  second  approach  is  a  much  more  original  treatment.  Indeed, 
there  are  no  results  in  the  literature  known  to  us  which  incorporate 
features  to  be  found  in  this  second  group  of  models.  In  this  approach, 
which  also  involves  dynamic  models,  there  is  assumed  a  priori  the  form 
of  the  probability  distribution  of  the  time  of  incidence  of  obsolescence 
and  the  parameter  values  for  this  distribution.  However,  the  parameter 
values  are  modified,  within  the  model,  by  Bayesian  procedures  involving 
a  functional  relationship  between  the  recent  history  of  demand  for  the 
line  item  and  the  parameter  values  of  the  obsolescence  distribution. 

This  basic  approach  is  adapted  in  several  ways  in  the  models  of 
this  sort  which  are  discussed  herein.  In  one  model  the  crucial  effect 
of  demand  history  appears  when  demand  is  at  a  level  of  zero.  In  another 
there  is  incorporated  a  Markov  process  which  deals  with  the  transition 
between  various  levels  of  demand.  In  a  third,  there  is  added  on  a 
feature  involving  the  dependency  of  demand  for  one  item  upon  the  de¬ 
mand  for  another  related  (complementary  or  substitute)  item. 

What  is  true  of  all  these  models  is  their  dynamic  character. 

They  are  true  dynamic  programs  which  optimize  with  respect  to  a  cost 
function  which  includes  all  cost  elements  of  an  inventory  system  ex¬ 
cept  relocation  of  existent  inventory  within  the  system,  and  repair  of 
items.  Thus,  this  family  of  models  does  not  deal  with  questions  of  how, 
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within  an  inventory  system,  the  inventory  should  be  dispersed.  Neither 
does  it  deal  with  repair  vs.  purchase,  nor  with  level-of-complexity  of 
items  (the  problem  referred  to  by  A.J.  Clark  as  the  "Inventory  Range 
Problem.") 

In  subsequent  chapters  these  models  will  be  dealt  with  in  de¬ 
tail.  In  addition,  there  will  be  some  discussion  of  the  results  of 
statistical  analysis  of  data  on  some  4000  aircraft  spares  which  have 
been  declared  obsolete  during  the  past  year.  Finally,  there  will  be 
some  discussion  of  explorations  which  were  undertaken  in  connection 
with  the  development  of  the  most  sophisticated  of  the  inventory  models 
which  use  the  properties  of  a  class  of  probability  distributions  that 
reproduce  under  the  operation  of  taking  the  conditional  distribution. 
The  use  of  elements  of  this  general  class  of  functions  in  these  models 
to  represent  the  probability  distribution  of  time  to  obsolescence 
leads  to  great  simplification  of  the  computing  problem.  This  is  so 
because  for  these  functions  the  parameters  are  computable  by  simple 
recursions  which  are  time- independent . 
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CHAPTER  III 


SURVEY  OF  INVENTORY  MODELS  WITH  EXPLICIT  TREATMENT  OF  OBSOLESCENCE 


Numerous  investigations  have  been  made  of  the  problems  of  inventory 
control  in  the  past  decade.  Although  most  of  them  have  dealt  with  the 
determination  of  optimum  policies  with  which  to  regulate  inventories  in 
the  absence  of  possible  obsolescence  of  the  inventories,  there  are  still 
a  number  of  studies  in  which  some  aspect  of  the  problem  of  obsolescence 
was  explicitly  treated.  This  chapter  is  essentially  an  annotated 
bibliography  of  the  studies  on  inventory  obsolescence.  The  bibliogra¬ 
phy  is  by  no  means  exhaustive,  but  every  study  which,  in  our  opinion, 
made  significant  contributions  to  this  area  of  inventory  control  study 
is  included. 

Results  which  emerged  from  these  inventory  obsolescence  studies 
are  varied  depending  on:  (1)  the  assumption  regarding  the  occurrence  of 
obsolescence;  (2)  the  nature  of  the  problem  that  the  investigator  has 
to  deal  with;  and  (3)  the  approach  adopted  for  obtaining  an  optimal 
solution.  In  order  to  bring  some  order  to  these  widely  disparate  re¬ 
sults  and  facilitate  further  discussion,  an  attempt  has  been  made  to 
arrange  the  material  to  be  discussed  in  a  systematic  way  in  which  the 
organization  is  based  on  the  problem- oriented  viewpoint  rather  than  on 
the  technique-oriented  viewpoint.  The  criteria  for  classification  are: 

1.  Deterministic  vs.  probabilistic  problem.  If  the  parameter 


which  specifies  the  occurrence  of  obsolescence,  such  as  the 
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length  of  time-to-obsolescence,  is  known  with  certainty, 
such  a  problem  is  referred  to  as  deterministic.  If  it  is 
known  in  terms  of  a  probability  distribution,  it  is  called 
the  probabilistic  problem.  A  problem  will  be  also  con¬ 
sidered  "deterministic"  if  obsolescence  enters  the  problem 
only  as  a  non- stochastic  function  of  any  inventory  parameter. 

2.  The  objectives  of  the  models.  Inventory  problems  can  be 

classified  according  to  the  objective  that  the  investigator 
sets  to  accomplish.  The  following  three  distinct  objectives 
are  noted  in  the  studies  to  be  reviewed  here: 

a.  The  main  objective  is  to  determine  the  optimum  economic 
lot  size  for  a  single  item  taking  cognizance  of  ob¬ 
solescence. 

b.  The  determination  of  the  optimum  inventory  for  a  single 
item  if  no  additional  procurements  are  to  be  made.  This 
optimum  inventory  is  commonly  referred  to  as  final  in¬ 
ventory. 

c.  To  determine  the  (s,S)  type  ordering  and  disposal  policies 
taking  account  of  obsolescence  cost. 

In  Table  I,  the  inventory  obsolescence  studies  are  classified  accord¬ 
ing  to  the  criteria  discussed  above  and  they  will  be  discussed  at  some 
length  in  the  order  indicated  in  each  entry  of  the  table. 


OBJECTIVE 


Table  I:  Classification  of  Inventory  Obsolescence  Studies 


Assumption  about  Occurrence  of  Obsolescence 


Deterministic 


Case  1 


Whitin  &3] 


Case  3 

Hadley  and  Whitin  £7] 


Case  5 


Barankin  00 
Fukuda  00 


Probabilistic 

Case  2 

Grass i  and  Gradwohl  [fj 

Case  4 

Hadley 

DO 

Mohan  and  Garg 

C11] 

S imps  on 

[12] 

Case  6 

Allen  and  Broida 

CO 

Barankin  and  Denny  [2] 

Ford 

CO 

Brown 

00 

|cQ| 

0«! 
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Case  1;  Deterministic  BOQ  Model 

The  solution  of  the  economic  order  quantity  problem,  according  to 
Whitin  [13],  has  a  rather  long  history.  It  is,  perhaps,  still  one  of 
the  most  frequently  used  inventory  formulae.  In  the  formula  for  de¬ 
termining  economic  order  quantities,  obsolescence  cost  is  treated  as  a 
constant  percentage  of  inventory  carrying  cost  which  increases  as  in¬ 
ventories  increase.  This  carrying  cost  is  to  be  balanced  against  those 
cost  items  such  as  quantity  discounts,  freight  differentials  which  de¬ 
crease  as  inventories  increase.  The  derivation  of  the  formula  is  as 
follows : 

Let 

Y  =  expected  demand  per  period 

Q  =  economic  order  quantity 

C  =  unit  inventory  carrying  cost  per  period.  (This  is 
made  up  of  material,  interest,  depreciation  and  ob¬ 
solescence  costs) 

S  =■  procurement  expense  per  order 

Then  total  variable  cost  (TVC)  per  period  is 

(3.1)  TVC  =  +  Q  S 

Upon  differentiating  (1)  with  respect  to  Q  and  setting  the  re¬ 
sulting  derivative  equal  zero,  we  obtain  the  minimizing  value  of  Q 
by  solving  the  derivative  for  Q. 


As  can  be  seen  in  the  above  formula,  if  the  risk  of  obsolescence 


should  increase,  this  will  be  reflected  in  a  higher  inventory  carrying 
cost  and,  consequently,  a  smaller  economic  order  quantity. 

In  addition  to  the  EOQ  model  described  above,  many  traditional 
inventory  models  treat  obsolescence  cost  as  a  part  of  carrying  cost. 
However,  when  obsolescence  becomes  a  major  cost  in  an  inventory  model, 
a  more  explicit  determination  of  obsolescence  costs  becomes  desirable. 
Case  2.  Probabilistic  EOQ  Model 

In  this  model  obsolescence  cost  is  formed  as  a  function  of  the 
obsolescence  rate  by  assuming  that  the  probability  of  obsolescence  at 
a  future  time  can  be  expressed  by  a  certain  distribution  function. 
Grassi  and  Gradwohl  £6^  have  obtained  a  probabilistic  EOQ  formula  by 
assuming  the  life  span  of  an  item  to  follow  an  exponential  density 
function,  f(t)  =  ne’^e,  where  \i  is  the  obsolescence  or  death  rate. 
Their  model  may  be  stated  as  follows: 

Let 

P  =  unit  order  cost 

D  =  the  sum  of  unit  material,  labor  and  overhead  costs 
S  =  setup  cost  per  order 
R  =  expected  demand  per  period 
B  =  the  safety  or  buffer  stock 
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I  -  unit  inventory  holding  cost  per  period 
P  *  the  inventory  charge  rate  per  period 
E  =  expected  unit  obsolescence  cost 
Then  the  sum  of  all  the  costs  associated  with  ordering  one  unit  of  in¬ 


ventories  is: 


C  =  P  +  I  +  E 


Each  component  of  the  total  unit  cost  C  is  given  by: 


P  =  D  +  ? 


P  may  be  regarded  as  the  value  of  one  unit  of  inventories. 


(3-4)  I  *  (B  +  |  )  (D  +  |  )  | 

B  +  ^  is  the  average  inventory  level  under  the  assumption  that  a 
quantity  will  be  ordered  when  the  inventory  level  reaches  the  safety 
level  B;  the  product  (B  +  7}  )  (D  +  ^  )  is  the  value  of  the  average 
inventory  stock.  The  inventory  charge  is  made  on  this  average  inventory 
stock  and  is  prorated  over  the  total  expected  requirement  for  one  period. 

The  expected  unit  obsolescence  cost  may  be  calculated  based  on  the 
following  consideration:  Suppose  a  lot  size  of  Q  is  ordered  for  an 
item  when  the  item  is  at  age  t  years  (unit  for  measuring  the  passage 
of  time  is  arbitrary).  The  lot  will  be  used  up  when  the  item  is  of  age 
T  +  R  year8,  It  is  then  of  interest  to  know  the  probability  of  obso¬ 
lescence  during  the  time  interval  ^  ,  given  non-obsolescence  prior  to 
t  so  that  the  expected  loss  due  to  obsolescence  during  the  time  inter- 
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val  of  length  ^  can  be  calculated. 

R 

Because  of  the  assumption  that  the  life  span  of  the  item  follows  an 
exponential  distribution,  the  probability  that  obsolescence  occurs  prior 
to  t  (or,  alternatively,  the  probability  that  the  item's  life  span  is, 
at  most,  as  long  as  t)  can  be  stated  as  follows: 


Pr (T  <  t)  =  F  (t)  =  l-e‘^ 

where  T  represents  the  age  at  obsolescence  of  that  particular  item  in 
question. 

The  conditional  probability  F  (T!  tQ)  of  obsolescence  in  an  addi¬ 
tional  time  T  after  experiencing  non-obsolescence  to  tQ  is  given 


G  (T) 


F(t 0+f)  -  F(to) 

~TTTT) 


Hence  the  conditional  density  function  is 


=  g  (f)  = 

The  expected  number  of  unit  lost  due  to  obsolescence  during  the  time 

Q 

interval  t  to  t  +  „  is: 

R 


=Jq  (B  +  Q  -  Rp  g (T) 


(3.5)  =  (B-^  )  (l-e‘^(Q/R))  +  Q 
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The  expected  unit  loss  due  to  obsolescence  becomes 


(3.6)  E  =  (D  +  S  )  (B-R/n)  (1- 

Q  Q 


-H(Q/R), 


The  optimum  economic  order  quantity  is  obtained  by  finding  the  mini¬ 
mizing  value  of  Q  in  the  following  expression 


P  +  I  +  E 


!>[«♦<.♦§>  I  -  ^ r"<Q/R>)] 


The  resulting  solution  given  by  Grassi  and  Gradwohl  is 


^RS  [l  +  |  (u  +  P)] 


(H  +  P) 


It  is  interesting  to  note  that  if  the  conditional  density  g(T)  follows 

a  rectangular  distribution,  i.e.,  obsolescence  is  equally  likely  to 

occur  at  any  moment  after  the  item  has  survived  up  to  t,  the  new  EOQ 

formula  is  exactly  the  same  as  (3.8)  except  the  second  term  in  the  de- 
Bu2 

nominator  vanishes. 

K 

In  using  the  formula  (3.8)  the  level  of  buffer  stock  (3  is  assumed 
to  be  determined  outside  of  the  model.  If  the  joint  determination  of  Q 
and  B  is  desired,  it  may  be  achieved  by  appending  penalty  cost  due  to 
stock  shortage  to  the  cost  function  (3.7)  and  minimizing  the  function 
with  respect  to  B  and  Q. 
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Case  3.  Final  Inventory  Model  with  a  Known  Obsolescence  Date 

The  problem  is  to  determine  the  optimal  amount  of  inventory  on  hand 
if  no  additional  procurements  are  allowed,  and,  furthermore,  a  date  of 
obsolescence  is  known  with  certainty.  This  inventory  is  commonly  known 
as  the  optimum  final  inventory.  The  main  feature  of  such  a  problem  is 
that  it  involves  two  types  of  cost;  one  is  incurred  at  a  fixed  point  in 
time  (obsolescence  cost)  and  another  is  a  function  of  time  (holding  and 
stockout  costs). 

Hadley  and  Whitin  £7^  formulated  a  model  which  deals  with  problems 
such  as  above  and  obtained  a  solution  assuming  the  demand  for  the  item 
to  be  Poisson  with  the  mean  demand  rate  independent  of  time 

Let  k  ■*  the  liquidation  loss  per  unit  at  the  time  of  obsolescence 

h  *  the  carrying  cost  per  unit  per  unit  time 

x  =  the  stockout  cost  per  unit  per  unit  time 

p(X;  Xt)  =*  the  Poisson  probability  distribution  that  the  demand  is 
exactly  X  units  in  a  time  period  of  length  t  with  X 
being  the  average  demand  rate 
T  *  time  of  obsolescence 
H  «  the  optimal  final  inventory 

The  expected  cost  of  holding  H  units  from  time  o  to  T  is  the  sum 
of  the  expected  costs  of  disposing  the  unused  amount  at  T;  carrying 
inventory,  and  stockouts  for  the  time  interval  o  to  T.  These  cost 
components  may  be  stated 
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The  expected  disposal  cost: 


H-l 

K  £  (H  -  X)  p(X  ;  X  T) 
X-0 


The  expected  carrying  cost: 

Since  the  carrying  cost  "rate"  at  time  t  is 

H-l 

h  £  (H  -  X)  p(X  :  X  t) 

X=0 

the  total  carrying  cost  for  the  interval  0  to  T  is 


(3.10) 


f  H-l 

‘  f  S 

J  X=0 


(H  -  X)  p(X  :  X  t)  dt 


the  expected  stockout  cost: 

Similarly,  the  total  stockout  cost  for  the  interval  0  to  T  is 


(3.11) 


£  (X  -  H)  p(X  :  X  T)  dt 
X=H 


The  expected  cost  of  holding  H  which  is  the  sum  of  (3.9),  (3.10) 
and  (3.11)  can  be  minimized  by  noting  the  identity 


P(X  }  \  t)  dt  =  -  £  p(u,  X  T) 

K  u  =  X  +  1 


This  model  is  neat  and  simple  to  apply.  However,  its  usefulness  may  be 
limited  for  the  following  reasons: 

1.  In  general,  it  is  not  advisable  to  assume  that  an  obsolescence 
date  can  be  specified  with  certainty. 
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2.  In  the  lifetime  of  an  item,  the  practice  of  repeated  procure¬ 
ments  may  be  more  common  than  the  situation  represented  by  the  model. 
Case  4:  Final  Inventory  Models  with  Probabilistic  Obsolescence  Date 


One  of  the  objections  to  the  model  in  Case  3,  is  in  regard  to  the 
assumption  of  known  obsolescence  date.  This  assumption  was  relaxed  in 
a  more  recent  model  developed  by  Hadley,  [8]  . 

He  assumed  that  the  date  of  obsolescence  can  be  described  by  a 
continuous  density  function  r(T).  Then  the  probability  that  the  item 
will  become  obsolete  in  the  interval  T  to  T  +  dT  is  r(T)dT.  With 
this  new  assumption  regarding  T,  the  expected  cost  of  holding  H  may 
be  modified  by  taking  the  expected  value  of  (3.9),  (3.10),  and  (3.11) 


with  respect  to  T. 


(3.12) 


E(H)  =  k  S  (H  -  X)  I  r (T)  p(X;  X.  T)  dT 
X=0  J 

T  0 

r  c  H‘1 

+  h  I  f  2  (H  -  X)  p(X;  \  t)  r  (T)  dT  dt 

J  J  X-0 


+  h  i  /  i 

00  T 

♦  «r  f  ; 

J  J  X=H 


(X  -  H)  p(X;  X  t)  r  (T)  dT  dt 


T=0  t-0 


Using  the  technique  described  in  the  earlier  paper  by  Hadley  and 
Whitin  [7],  the  cost  function  (3.12)  can  be  minimized.  Hadley  [fe]  has 
also  developed  another  model  in  which  only  a  finite  number  of  possible 
dates  of  obsolescence  is  assumed  and  the  probability  of  each  is  specified. 

Another  type  of  final  inventory  model  was  developed  by  Simpson  [12] 
and  Mohan  and  Garg  [11].  Their  model  is  distinguished  from  Hadley's 
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model  in  that  only  a  knowledge  of  the  average  annual  demand  is  assumed, 
and  it  will  first  calculate  the  optimal  economic  retention  period  from 
consideration  of  the  balancing  of  carrying  charges  and  the  cost  of 
disposal.  More  specifically,  the  cost  of  retaining  a  unit  value  of  the 
stock  is  weighed  against  the  cost  of  disposing  one  unit  value  of  the 
stock  now  and  procure  it  again  at  some  later  stage.  It  is  clear  that 
units  will  be  added  to  the  retention  stock  as  long  as  the  latter  cost 
exceeds  the  former.  On  the  other  hand,  the  retention  stock  will  be 
reduced  if  the  former  costs  exceed  the  latter.  The  equilibrium  is 
reached  when  the  cost  of  storing  the  marginal  unit  is  exactly  equal  to 
the  cost  incurred  for  not  storing  that  unit. 

Their  model  may  be  stated  as  follows: 

Let  U  =  average  annual  demand 

x  =  the  number  of  years  for  which  a  stock,  say  N,  will  meet 
average  demand.  Note  the  relation  xU=N 
D  =  fraction  of  unit  value  of  material  which  will  be  realized 
in  disposal  sales 
i  =  interest  rate 

r  *  annual  storage  cost  of  material,  expressed  as  a  fraction 
of  the  unit  value  of  material 


cQl 


iQ 
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CD  «*  total  cost  of  retaining  the  unit  value  of  stock 
R 

Cp  *  total  cost  of  disposing  the  unit  value  of  stock 
Ft  =  the  probability  that  the  item  becomes  obsolete  prior  to 
the  (t+l)th  year 

The  retention  cost  CD  consists  of  two  types  of  costs,  (1)  the 

R 

storage  cost  for  x  years,  (2)  the  obsolescence  cost. 

The  first  type  of  cost  may  be  stated  as 

x  x-t+1 

r  £  (1  -  F  )  (1  +  i)x 

t=l 

(1  -  F  )  is  the  probability  of  non-obsolescence  before  the  (t+l)th 

x- 1 

year,  and  (1  +  i)  is  the  compound  interest  charge  for  storing 

the  unit  value  of  stock  for  x-t+1  years. 

The  obsolescence  cost  is  the  product  of  the  unit  value  of  stock  and 
the  probability  of  obsolescence  Fx- 
Hence 

(3.13)  C  =  F  +  r  £  (1  -  f  )  (1  +  i)x't+1 

R  X  t=l 

For  each  unit  value  of  stock  disposed  now,  the  cost  equivalent  to  the 
unit  value  will  be  incurred  in  order  to  procure  it  again  after  x  years 
assuming  there  is  no  price  change  in  the  future.  However,  for  each  unit 
value  of  stock  disposed,  a  salvage  value  of  D  is  acquired.  This  amount 
should  increase  at  compound  interests  for  x  years  and  should  be  con¬ 
sidered  as  a  credit.  Therefore,  the  disposal  cost  is 

CD  «  1  -  D(1  +  i)X 


(3.14) 
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The  optimum  retention  period  x  is  found  by  solving  the  equation 
CR  "  CD  =  °* 

The  only  difference  between  the  result  obtained  by  Simpson  and  that 
obtained  by  Mohan,  et  al  is  that  the  former  assumes  that  the  probability 
of  obsolescence  follows  a  uniform  distribution  while  the  latter  uses  a 
normal  distribution. 

Simpson  gives  a  detailed  account  of  the  advantages  and  the  short¬ 
comings  of  this  retention  and  disposal  formula  from  the  practical  point 
of  view.  His  conclusion  is  that  despite  the  fact  that  many  rigid 
assumptions  are  required  to  derive  this  formula,  it  provides  a  workable, 
practical  decision  rule. 

Case  5:  Optimal  Ordering  and  Disposal  Policies  with  Known  Obsolescence 
Date 

When  the  items  in  inventory  are  known  to  become  obsolete  in  the 
future,  it  is  advisable  to  consider  the  disposal  of  surplus  items  as 
well  as  the  procurement  of  stock  to  meet  the  future  demand. 

Consider  the  following  situation  in  which  an  item  will  no  longer  be 
used  after  a  certain  date  and  the  stock  on  hand  on  that  date  will  have  to 
be  disposed.  Meanwhile,  prior  to  the  occurrence  of  obsolescence,  the 
stock  level  of  this  item  is  reviewed  periodically,  say,  at  the  beginning 
of  a  finite  number  of  equal  time  intervals;  and  one  of  the  following  de¬ 
cisions  is  made 

1.  Procurement  of  additional  stock 

2.  Disposal  of  excess  stock 

3.  No  procurement  and  no  disposal 
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For  such  a  multi-stage  decision  process,  the  technique  of  dynamic 
programming  has  been  known  to  furnish  a  suitable  framework  for  analysis. 
In  the  following  the  recursion  relation  of  a  dynamic  programming  model 
appropriate  to  the  situation  described  above  is  stated. 

Let  us  introduce  the  following  definitions: 
s  =  demand  per  period 

x  =  "initial"  inventory  level;  an  inventory  level  at  the  beginning 
of  a  period  before  the  order  or  disposal  is  made 
y  =  "starting"  inventory  level;  an  inventory  level  at  the  beginning 
of  a  period  immediately  after  the  order  or  disposal  is  made 
r(z)  =  the  cost  of  ordering  z  units 

d(z)  =  the  cost  of  disposing  z  units.  If  the  disposal  item  has  a 
salvage  value,  this  cost  assumes  a  negative  value. 
h(z)  =  the  cost  of  holding  z  units  of  inventory  for  one  period 
p(z)  =  the  cost  per  period  of  having  demand  z  units  greater  than 
inventory 

p  =  a  discount  factor 

It  is  assumed  that  there  is  no  delay  either  in  the  order  or  in  dis¬ 
posing  of  stock.  Suppose  there  are  n  periods  before  obsolescence 
occurs.  The  policy  to  be  considered  is  of  the  following  simple  forms: 
for  each  period  a  starting  inventory  y  is  specified  for  each  value  of 
initial  inventory  x.  Associated  with  each  y  is  a  certain  expected 
cost.  The  problem  is  to  minimize  this  expected  cost  by  suitable  choice 
of  y. 


IcQ 

□  « 
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Let  us  define 


(x)  ■  expected  total  cost  for  a  process  which  has  k  periods 
remaining  and  starting  the  process  with  x  units  of 
initial  inventory  with  an  optimal  ordering  and  disposal 
policy. 

Since  this  total  expected  cost  is  composed  of:  (1)  a  cost  of  using 
an  optimal  policy  in  this  period,  (2)  a  discounted  future  cost,  we  have 
the  following  recursion  relation 

(3.15)  CR(x)  -  Min  ^r(y-x)  +  d(x-y)  +  E[h(y-s)  +  p(s-y)J  +  p  E  ck_1(y*8)}' 

y>0  for  k  «  2,3,...,n. 

where  the  ordering  cost  and  disposal  cost  functions  are  zero  for 
negative  arguments.  E  denotes  expectation  and  the  expectation  opera¬ 
tion  is  with  respect  to  demand  a. 

When  there  is  only  one  period  remaining  in  the  process,  the  total 
expected  cost  is  simplified  to  the  following: 

(3.16)  Cj  (x)  -  Min^r(y-x)  +  d(x-y)  +  E  £h(y-s)  +  p  (a-y)J 


Together,  (3.15)  and  (3.16)  enable  us  to  determine,  successively 
for  k“l,2, . . . ,n,  the  optimal  value  of  y  for  each  x  and  the  corre¬ 
sponding  minimum  expected  cost  C^(x). 

This  is  the  type  of  model  formulated  by  Barankin  and  Denny  ^2j  and 
Fukuda  [5] .  Although  the  cost  functions  as  well  as  the  demand  per  period 
are  considered  to  be  independent  of  time  in  the  above  model,  the  same 
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technique  still  applies  even  if  they  should  become  time  dependent.  The 
assumption  of  no  time  lag  in  delivery  can  be  also  relaxed  and  the  same 
technique  can  be  used  as  long  as  a  time  lag  is  known  with  certainty.^ 

Case  6;  Optimal  Ordering  and  Disposal  Policies  with  Probabilistic 
Obsolescence  Date 

Barankin  and  Denny  [2]  and  Ford  £4]  have  extended  the  model  de¬ 
scribed  in  Case  5  to  the  situation  in  which  the  information  about  obso¬ 
lescence  is  available  in  the  form  of  a  probability  statement,  e.g.,  the 
probability  of  obsolescence  in  a  given  period,  say  the  kth  period,  is 
d^.  Presumably,  Sd^  =  1.  Then  the  conditional  probability  a^  that 
obsolescence  occurs  in  period  k,  given  non-obsolescence  in  periods 
n,n-l,4..,k  +  1,  may  be  calculated  as  follows: 


1  d- 

1=1  1 

With  a  knowledge  of  these  conditional  probabilities  for  k=l,2,...,n 
only  a  simple  modification  of  the  recursion  relation  (3.15)  is  needed  to 
give  us  a  desired  inventory  model. 

Ck(x)  =  Min^c(y-x)  +  E  P  (s-y)  +  E  h(y-s)  +  akd(x-y)+(l-ak>  p  E  Ck+1(y-s) 
yX)  for  k  =  2,3,...,n 


See  H.  Scarf,  "(SjS)  Policy  in  Dynamic  Inventory  Problem, "Mathematical 
Methods  in  the  Social  Science,  Stanford  University  Press,  Standord, 
California. 
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Since  the  conditional  probability  of  obsolescence  in  the  last  period  Sj 
is  equal  to  unity,  tne  total  expected  cost  function  for  k-1,  is  the  same 
as  (3.16).  Again,  starting  with  k-1,  we  can  successively  calculate  a 
set  of  optimal  policies  and  corresponding  expected  total  costs  for  each 
period. 

A  more  detailed  discussion  regarding  the  modeling  of  this  type  of 
inventory  problem  as  well  as  some  characteristics  of  optimum  policies 
is  presented  in  Chapter  IV  of  this  report. 

In  all  the  dynamic  inventory  models  discussed  so  far,  recursion  re¬ 
lation  such  as  (3.15)  is  always  characterized  completely  by  two  state 
variables,  x,  the  stock  on  hand,  and  n,  the  number  of  remaining 
periods.  Brown  considered  a  more  realistic  situation  in  which  the 
latter  state  describing  variable  is  given  in  terms  of  a  probability  dis¬ 
tribution.  This  introduction  of  a  stochastic  variable  for  describing  the 
nature  of  process  does  not  unduly  complicate  the  calculations.  Results 
of  this  type  of  dynamic  inventory  model  can  be  found  in  Chapter  VI  of 
this  report. 

In  the  study  made  by  Stanford  Research  Institute  for  BuSandA, 

Allen  and  Broida  [l]  considered  the  problem  of  minimizing  the  unit-weeks 
of  system-wide  shortages  subject  to  the  total  variable  procurement 
budget  (converted  into  the  shortage  equivalent)  by  suitable  choice  of 
non-negative  order  quantities.  Since  the  stock  level  of  each  item  in 
the  system  is  reviewed  every  quarter  and  an  order  is  placed  if  necessary 
to  raise  the  stock  level  to  a  certain  level  which  is  most  desirable  from 
the  standpoint  of  minimizing  the  shortage  risk,  Allen  and  Broida  set  out 
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to  determine  what  should  be  considered  the  most  desirable  level  with  a 
fixed  amount  of  resources  available  to  the  system. 

The  suggested  procurement  rule  is  as  follows: 

Let  R  =  the  most  desirable  inventory  level  (determined  by 
the  model) 

I  =  the  sum  of  stock  on  order  and  on  hand 

(1)  If  R  -  I  <  0,  no  order  is  placed 

(2)  If  R  -  I  >  0,  order  the  difference  (R-I)  or  the  economic 
order  quantity- -whichever  is  greater 

This  procedure  may  be  considered  a  variant  of  the  usual  (s,S)  policy. 
The  most  interesting  feature  of  their  study  is  that  a  budget  constraint 
was  explicitly  taken  into  account  in  deriving  an  optimal  inventory  regu¬ 
lating  procedure. 

Since  the  derivation  of  their  procurement  rule  is  rather  lengthy, 
it  will  not  be  reproduced  here  except  to  indicate  how  obsolescence  costs 
are  incorporated  into  their  model. 

They  assumed  that  the  probability  of  obsolescence  at  a  future  time 
to  be  expressed  by  an  exponential  distribution 

I-'* 

vihere  Y  i-8  the  obsolescence  factor  for  the  item  in  question,  the 
probability  of  non-obsolescence  is  then  e  .  The  probability  of  non¬ 
obsolescence  at  time  t  is  multiplied  by  the  expected  number  of  units 
short  at  that  time.  This  yields  the  expected  number  of  units  short  at  time 
t  when  no  obsolescence  occurs.  From  this  last  quantity,  the 
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expected  number  of  unit-weeks  of  shortage  is  obtained  by  integrating  out 
t.  Similarly,  in  estimating  the  budget  expenditure  needed  for  placing 
an  order  of  a  certain  size,  it  is  weighed  by  the  probability  of  non¬ 
obsolescence* 

When  obsolescence  is  treated  in  the  manner  described  above,  it  has 
a  net  effect  of  reducing  the  amount  of  budget  expenditure  available  for 
procurement  purposes;  at  the  same  time  it  lowers  the  risk  of  shortage 
because  when  the  item  has  already  become  obsolete  there  is  no  question 
of  shortage.  It  is,  therefore,  difficult  to  generalize  the  over-all 
effect  of  the  above  modeling  of  obsolescence  on  the  resulting  procure¬ 
ment  policy  without  examining  in  more  detail  a  specific  parameter  value 
used  to  describe  the  occurrence  of  obsolescence  for  a  problem  in  question. 


0=1 
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CHAPTER  IV:  A  DYNAMIC  PROCUREMENT-DISPOSAL  MODEL 

A.  Introduction 

In  this  chapter  is  developed  a  dynamic  procurement-disposal  in¬ 
ventory  model  which  assumes  that  the  number  of  time  periods  until 
obsolescence  and  the  demand  for  each  line  item  are  random  variables 
with  known  probability  distributions.  While  the  model  assumes  that  the 
demand  distribution  is  known  for  an  item  for  each  time  period  of  its 
lifetime,  it  allows  this  distribution  to  vary  from  period  to  period. 
Thus,  diminishing  demand  toward  the  end  of  the  item's  lifetime  may  be 
reflected.  In  addition,  although  the  probability  distribution  of  time 
to  obsolescence  is  assumed  known,  it  too  can  vary  from  period  to 
period.  Also  assumed  known  are:  holding  cost  per  unit  per  time  period; 
stockout  cost  per  unit  per  time  period;  unit  price;  fixed  reorder  cost; 
and  disposal  cost  per  unit. 

Each  of  the  known  costs  referred  to  above  may  be  varied  from  peri¬ 
od  to  period.  Thus,  any  considerations  which  suggest  that  costs  will 
change  in  some  fashion  in  the  future  can  be  taken  into  account  in  com¬ 
puting  policy. 

A  further  restriction  is  imposed  solely  for  reasons  of  computing 
simplicity.  This  is  that  lead  time  is  equal  to  one  inventory  period. 

The  program  is  capable  of  modification  to  allow  lead  time  to  be  a  mul¬ 
tiple  of  the  inventory  period,  but  for  purposes  of  numerical  analysis  it 
was  felt  that  at  this  time  the  great  increase  in  computing  would  not  be 


justified 
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B •  Some  General  Remarks 

It  should  be  noted  that  the  dynamic  program  operates  backwards 
in  time.  It  begins  by  determining  an  optimal  policy  for  some  "last"'*' 
time  period,  proceeds  from  that  to  an  optimal  policy  for  the  next  to 
last  period,  and  so  on.  Thus,  for  an  obsolescence  date  of  twenty 
periods  into  the  future,  the  twentieth  from  the  "last"  period  is  the 
initial  period  and  the  best  policy  for  that  period  is  the  best  initial 
policy.  A  "best  policy"  for  a  period  is  defined  by  the  program  as 
follows.-  it  is  a  policy  which  specifies  for  each  feasible  stock  level 
left  on  hand  at  the  end  of  the  previous  period  what  the  level  ought  to 
be  brought  to  at  the  beginning  of  this  period  in  order  that  the  sum  of 
all  costs  from  this  period  on  shall  be  minimal. 

The  program  is  based  on  the  fact  that  at  the  start  of  any  time 
period  the  optimal  stock  level  depends  only  on  the  inventory  left  from 
the  previous  period  and  not  on  how  that  inventory  level  came  about; 
that  is,  the  past  influences  the  present  through  a  single  number,  the 
inventory  level  at  the  end  of  the  immediately  preceding  period.  For 
that  reason,  the  program  is  able  to  specify  the  future  cost  associated 
with  starting  a  time  period  with  each  possible  inventory  level,  without 
reference  to  what  inventory  will,  in  fact,  be  left  over  from  the  pre¬ 
ceding  period  or  what  policy  was  used  in  the  earlier  periods.  It  only 


This  "last"  period  is  determined  by  examining  the  probability  distri¬ 
bution  of  obsolescence  occurring  at  time  t (t=0, . . . , f , . . .,  ») .  We 
find  a  T  such  that  the  probability  of  obsolescence  occurring  after  T  is 
less  than  some  arbitrary  e.  T  is  that  "last"  period. 


-35- 


needs  to  assume  that  an  optimal  policy  is  used  in  all  later  periods. 
This  fact,  plus  a  knowledge  of  what  policy  is  optimal  during  the  first 
period  after  obsolescence  (specifically,  to  dispose  of  all  stock  on 
hand)  provides  the  recursive  procedure  of  the  dynamic  program. 

The  program  produces,  then,  a  policy  for  each  time  period.  Each 
policy  is  used  to  calculate  a  preceding  policy,  and  one  entry  of  each 
(the  best  starting  stock  when  the  inventory  left  from  the  previous 
period  is  zero)  is  also  used  to  provide  an  optimal  initial  order  for 
each  possible  number  of  periods  to  obsolescence. 

Furthermore,  since  the  full  set  of  order  policies  include  optimal 
stock  levels  when  there  is  already  some  stock  left  over  from  the  past, 
the  model  is  fully  capable  of  suggesting  reorder  decisions  as  well  as 
initial  order  amounts. 

C .  A  Specific  Formulation  of  the  Dynamic  Program 
1.  Definitions 

a.  Cost  of  changing  stock  level 

Let  ja,bi  be  a  known  function  denoting  the  cost  in 

period  i  of  converting  a  stock  level  of  "a"  units  into 
a  stock  level  of  "b"  units.  (a,b  >  0). 

It  is  postulated  that  y^  has  the  following  properties: 

(4.1)  Y  ja,b.  +  Y I h,c  |  >  Yfa,cl 

(4.2)  Y;a+k,  b+kj  =  Y'a;hj  for  all  k  for  which  >  0 

(4.3)  Y  ja,  b |  >  0,  y\a,a  =  0 
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The  purchase  price  of  an  item  is  not  included  in 
reorder  cost  (V|0>+|);  it  is  contained  instead  in 
disposal  cost  (y|o,-j). 

Making  use  of  (4.2)  we  may  define  y|a,bj  for  a,b 
including  negative  values,  and  may  then  define  0^(a) 

=  y|o,aj  =  yj\>,  bfaj  for  all  a  =  0,  +  1,  +  2,  .... 
0^(a)  shall  be  used  as  the  input  quantity. 

b.  Shortage  cost 
Let 

Cost1  (r-k)  =  Cost  in  ith  period  associated  with 
out  — 

an  incoming  order  r  which  exceeds  the  beginning 
inventory  k,  r-k  =  1,  2,  3  .... 

This  cost  will  be  specified  as  an  analytic  function 
rather  than  a  table  of  values  in  the  computer  program. 

c.  Storage  cost 

Cost1  (r,k)  =  Storage  costs  in  ith  period  associated 
with  an  initial  stock  of  k  and  an  incoming  order  of 
r.  If  this  cost  is  to  be  charged  at  the  end  of  a  peri¬ 
od,  it  can  be  regarded  as  a  function  of  k-r,  the  ex¬ 
cess  of  inventory  over  demand.  Since  this  cost  is 
defined  for  each  time  period,  any  discounting  of  the 
future  or  interest  charges  on  investment  will  be  in¬ 


cluded  in  it. 
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d.  Issue  cost 
Let 

Cost^sg(r)  =  cost  of  issuing  r  units  in  the  ith 
period. 


r  =  1,  2,  3,  .... 

Note  that  this  cost  must  include  the  unit  purchase 
cost  in  order  for  the  program  to  operate  properly, 
e.  Probability  Distribution  of  Demand 

Let  Pr  (X^  =  r)  =  Probability  that  the  demand  in 
period  i  is  exactly  r  units. 


r  =  0,  1,2,  .... 

Probability  Distribution  of  Obsolescence 
Let  a^  =  Conditional  probability  that  obsolescence 
occurs  in  period  i,  given  that  it  has  not  occurred 
in  periods  N,N-1,  ...,  i+1. 

More  specifically 
d, 


ai  =  T 


Z  d 

j  =  l 


j 


where  d^  is  the  probability  that  obsolescence  occurs 

n 

in  period  i,  and  £  d.  =  1. 

j=l  1 

Expected  future  costs 
Let 

jeT(K)=  Expected  cost  from  i  on,  if  the  initial  stock 
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level  in  the  ith  period  is  K:  K  =  0,  1,  2,  . ..,  and 
an  optimal  policy  is  used  in  period  i  as  well  as  in 
each  of  the  future  periods. 

2.  Induction 

Since  we  have  an  optimal  policy  for  time  period  zero  (i.e., 
if  there  are  K  units  in  inventory,  dispose  of  K  units) 
and  a  set  of  expected  future  costs  for  period  zero, 

0O(K)  =  eQ(-K)  for  K  =  0,  1,  2,  .... 

the  entire  dynamic  program  can  be  specified  by  describing  a 
procedure  for  calculating  policy  and  future  cost  in  period  i 
from  policy  and  future  cost  in  period  i-1.  Note  that  period 
i  is  the  time  period  which  is  i  periods  earlier  than  the 
date  of  obsolescence. 

Assume  now  that  0^  ^ (K)  has  been  calculated. 

(4.4)  Then: 

MK)  =  2  Pr(x  =r)  •  t.(K,r) 

i  r=0  1 

And  f^K,  r)  is  given  by: 

(4.5)  t.(K,r)  =  Cost  ^  (r-K)  +  Cost  gtQ  (K, r)  +  Cost  *sg  (r) 

+  (1-a  )  M  (r-K)  +  fl^r-K)  if  r  >  K 

=Cost  \  (K,r)+Cost  ]  (r)+(l-a  )M(r-K)+a  9  (r-K) 

SCO  1SS  1  i  i 

if  r  <  K 


cQ 

H« 
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H(r-K)  is  calculated  as  follows 
For  given  h,  h  =  0,  +  1,  +  2,  ... 

(4.6)  Define 

M(h)  =  min  1"  0 .  (k)  +  0.(h+k)l 

k>0  L  1  1  ] 

and  define  S^(h)  =  the  value  of  k  at  which  M(h)  is 
achie ved 

S^(h)  is  interpreted  as  the  ith  period  policy,  i.e., 
if  the  initial  stock  for  period  i  is  K  and  =  r, 
then  S^(r-K)  is  the  amount  with  which  to  start  period 
(i-1). 

Note  that,  in  (4.5),  the  expression  for  ^  involves 
all  the  costs  actually  incurred  in  period  i,  plus  the 
expected  cost  associated  with  the  ending  inventory, 

M(r,K),  if  the  process  is  to  continue;  and  the  disposal 
cost,  0(r,K),  if  the  process  is  to  terminate.  The  last 
two  costs  are  weighted  according  to  the  probabilities  of 
the  events  which  will  generate  these  costs. 

Some  remarks  regarding  the  calculation  of  M(h)  are 
in  order.  Assume  we  start  with  K  >  0,  there  is  an  order 
of  r  >  0  and  we  end  the  period  with  h  >  0  in  stock. 
Aside  from  costs  associated  with  K  and  r  above  which 
already  occur  explicitly  in  i|/\, 


we  also  have 
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CASE  I  r  >  K 

Additional  cost  =  ^ (k)  +  yi[K-r,kj=  0±  1  (k) 

+  9^  (k+r-K)  and  to  minimize  this  is  to  reorder  so  as 
to  bring  stock  level  to  S^r-K)  units,  and  additional 
cost  has  become  M(r-K) . 

CASE  II  r  <  K 

Additional  cost  =  0^  ^(k)  +  Yi[k-r,k  and  same 
minimization  as  in  Case  I. 

Note  that  the  argument  of  M(h)  is  non- positive. 

D.  Some  Remarks  on  Optimal  Policies  of  Dynamic  Procurement-Disposal 

Model 

A  number  of  computer  runs  were  made  in  order  to  study  the  effect 
of  parameter  variations  on  the  optimal  procurement  and  disposal  poli¬ 
cies.  Detailed  results  of  these  runs  are  presented  in  Appendix  A. 

Within  the  range  of  parameter  variations  explored,  it  was  found 
that  when  a  set-up  cost  is  involved  in  ordering  and  disposing  of  the 
stock,  the  optimal  ordering  and  disposal  policies  are  of  an  (s,S) 
type.  That  is  to  say  there  are  four  uniquely  determined  integers 
Kq  <  K-^  <  K2  <  such  that  if  the  initial  inventory  level,  say  H, 

is  less  than  K^,  an  order  is  placed  to  bring  the  inventory  level 
up  to  K^;  if  the  inventory  condition  is  Kq  <  H  <  K^,  no  action 
is  taken;  if  the  inventory  level  exceeds  K^,  it  is  disposed  down  to 
K^.  It  was  observed  that  these  policy  parameters  (K's)  are  related 
to  the  expected  total  inventory  cost  function  0(K)  (See  pp.  )  in 
the  following  manner: 

_  s®i 


q« 


-41- 


u  ”“o  1  2  3 

Fig.  1:  A  Typical  Expected  Total  Inventory  Cost  Function 


K,  ’inventory 
3  Level 


a  =  fixed  ordering  cost 

b  •  salvage  value  obtained  from  disposing  one  unit  of  the  stock 
c  *  unit  storage  cost 
d  »  fixed  disposal  cost 

1)  is  an  integer  value  of  K  which  minimises  0(K) 

2)  Kq  is  the  largest  integer  such  that 

0<KO)  >  0(Kl)  +  a  and  Kq  <  Kx 

3)  Kj  is  an  integer  such  that 

0'  (Kj)  <  b<0^(K1)  and  K2  > 

|  I 

0,  (K^)  and  0+  (Rj)  denote  the  left-hand  and 
right-hand  derivatives  at  respectively. 
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4)  is  an  integer  such  that 

C(K3-K2)  >  d  and  K3  >  K2 

These  findings  seem  to  make  sense  for  the  following  reasons: 

1)  Since  the  objective  of  the  model  is  to  minimize  the 

total  expected  inventory  cost,  it  is  natural  to  bring 
the  inventory  level  up  to  K^. 

2)  If  a  set-up  cost  is  involved  in  ordering,  the  cost  func¬ 
tion  in  effect  between  KQ  and  is  the  dotted 

line  above  the  0(K)  function  in  Fig.  1.  It  is  then 
most  economical  to  order  from  K^. 

3)  Suppose  there  is  no  fixed  disposal  cost.  As  the  initial 
inventory  level  increases,  the  total  expected  inventory 
cost  will  also  rise  because  of  an  increase  in  shortage 
cost.  At  some  point,  namely  K2,  it  becomes  cheaper 

to  dispose  an  additional  unit  of  the  stock  than  to  re¬ 
tain  it  since  in  this  model  a  salvage  value  is  attached 
to  a  unit  disposed. 

4)  If  a  set-up  cost  is  incurred  for  disposing  one  unit  of 
the  stock,  then  the  cumulative  storage  costs  must  exceed 
the  fixed  disposal  costs  before  a  disposal  becomes 
economically  attractive. 

Although  an  (s,S)  type  policy  seems  to  characterize  all  the  op¬ 
timal  policies  obtained  in  this  study  (see  Appendix  A),  even  when  a 
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constant  fixed  cost  is  involved  in  ordering  and  disposing  of  the  stock, 
ii  should  be  noted  that  such  a  policy  cannot  be  optimal  in  all  cases. 
However,  it  is  probably  safe  to  conclude  that  in  most  practical  situa¬ 
tions  an  (s,S)  policy  will  produce  a  near  optimal  solution. 


jcQ 

Oh 
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If,  in  fact,  Che  item  then  survives  until  t^,  the  probability  that  it 
will  survive  until  t 2,  t^  <  t^  <  T  would  be 


F2(t2)“Fl(t2ltl)’ 


F„<t  )-F  ,  (t  1 1  )  - 

n  n  n~ 1  n*  n-1 


r2 

1  f  dx 

*1  *2  -  fcl 

{H~rv 

*  j  * 

•0  0** 


tn~tn-l 

T‘t0'tn-1 


Clearly,  this  recursion  is  time-dependent.  Similarly,  if  we  assume 
that  log  t  is  Cauchy-distributed,  we  find  that 

.  P tan  1 (log  t  -H)  -  tan*1  (log  t  .-U)  "I 

F  (t  |t  )  -  i  - 2 - Szl - 

n  n1  n-1  « I  _ ,  ,  I 

Ll+2«(1+|)  -  tan  (log  tn*n)-tan  (log  tnl-H)J 

equally,  a  time-dependent  recursion. 

We  now  consider  some  examples  of  distributions  which  are  reprodu¬ 
cible  under  the  operation  of  taking  the  conditional.  The  simplest 

example  of  this  property  is,  of  course,  furnished  by  the  exponential 

-fct 

distribution,  with  probability  of  survival  to  age  t  equal  to  e  ; 
after  survival  to  age  t^  to  the  probability  of  surviving  for  a 
further  interval  of  t  is  again  e  .  We  present  three  related 
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farrili.fr  s  of  mortality  distributions  (including  the  exponential  as  a 
special  case)  which  exhibit  reproducible  conditional  distributions. 

If  the  cumulative  distribution  function  F(t)  is  represented 
in  terms  of  the  age- specific  death  rate  X(t),  we  have 


F  (t )  =  1  -  exp 


-  t 

-  J  \(u) 

-  0 


du 


The  conditional  probability  F(  *'  1 1^)  of  surviving  an  additional 
time  T,  after  surviving  to  tQ,  is  given  by 


F(t0+T)  -  F(t0) 

5^) 


=  1  -  exp 


tQ+T 

-  J  X(u) 


du 


I^-fc  pa'ticular  families  described  below  result  from  assuming  that 
F(T|tq)  is  related  to  F(t)  by  a  transformation  of  the  form  t— ^kT 


and 


W  X 


If 


is  differentiable  no  other  families  can  be 


reproducible  ander  the  same  transformations  (demonstration  omitted) 
Case  I:  ).(t)  =  ae6t,  k  =  1,  fc  =  e^0 


1  -  exp 


=  1  -  exp 


-j  X(u) 


du 


=  1  -  exp 


[■“/ 


*o+T 

***** 


a/  P(VT>  **0 


>] 
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=  1  -  exp 


=  1  -  exp 


(*o 


Pt°  f  PU  A 

-O.e  I  e.  du 
0 


Case  II 


\(t)  =  a(l+pt)Y;  k  =  — i—  ,  h  =  (1+pt n)Yfl 


1  -  exp 


t0+  ^ 

■  J  X(u) 

- 


du  I  =  1  -  exp 


(l+£u) r  du 


Ha:  Y  4=  -1 


1  *  exp 


[  r 

-a  J  (H 


(l+P'j)  y  du 


=  1  '  exp  pS=I)  l(l  +  p  (to+r))Yfl  '  (1+pto)Vfl 


=  1  ■  £Xp  '  p<?ni  afpt0)Vfl  {(1  +  )Vfl  _1} 


(5.2)  =  1-exp  -a(l+ptQ) 


i+pio 

i*+1J  (l+pu)Y  du 
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Ilb:  y  =  -1 


1  -  exp 


r0+T 

-a  J  (1 

-  tn 


(1+pu)  du 


a  ,  1+P<to+r> 
" 1 '  exp  '  p  ln  nsptj — 


1  -  exp  p  In  (  1  +  l+ptQ  ^ 


T 

1+Pfc0 


*.5  3)  =  1-exp 


^PL0 

-a  J  (1 


+  pu)  du 


Cumulative  distribution  functions 


Integrating  out  the  age-specific  death  rate  in  (5.1),  (5.2), 
and  (5.3),  and  reparameterizing  slightly  yields  the  following  cumu¬ 
lative  distributions: 


Case  I:  F(t)  =  1-e 


,  bt  .. 
-a(e  -1) 


The  conditional  reproducibility  of  this  function  may  be  demon¬ 
strated  as  follows: 

Let  F (t  ^ 1 1 q)  denote  the  probability  of  surviving  to  t^  given 
that  it  has  survived  up  to  t^. 
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bt. 


bt. 


M>-  1 +  .-‘O"0-1) 

btr, 


e-a(e  -1) 
bt0  . 


=  1-e 


By  defining  a  new  variable 


-ae  (e 


1) 


T=  C1  *  C0 


and  new  parameters 


bt 


0 


a1  =  ae  and  b’  =  b, 

We  obtain  a  new  distribution 

G(T)  -  l-e-'<b"r-l> 

which  has  the  identical  form  as  the  original  one. 

Similarly,  the  reproducibility  of  the  distributions  in  Cases  Ha 
and  lib  can  be  demonstrated. 


Case  Ila:  F(t) 


, [o+b=)c-i] . 


a'  =  a(l+bt0)  , 
b 


b'  = 


1+bt, 


Case  lib:  F(t)  =  l-(l+bt)  ;  b' 


1+bt0  ’ 


Note  that  Cases  I  and  lib  are  both  limiting  cases  of  Ila.  The  expon¬ 
ential  distribution  is  also  a  limiting  case  of  Ila. 
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= 


bt,  bt„ 

l_e-a(e  1‘1).  i  +  e‘a(e  °*1) 
bt0 

-a(e  -1) 


V  b(tfto) 


=  1-e 


-ae  (e 


By  defining  a  new  variable 


and  new  parameters 


r=ti-to 


a'  =  ae  and  b'  =  b, 


We  obtain  a  new  distribution 


G(T) 


=  i.e-a'Cb'T-l) 


which  has  the  identical  form  as  the  original  one. 

Similarly,  the  reproducibility  of  the  distributions  in  Cases  Ila 
and  lib  can  be  demonstrated. 


Case  Ila:  F(t) 


=  l-e-a  [d+bt:)c-l]  . 


a'  =  a(l+bt0)  , 


Case  lib:  F(t)  =  l-(l+bt)'c;  b'  =  tt 


1+bt0  ’ 


Note  that  Cases  I  and  lib  are  both  limiting  cases  of  Ila.  The  expon¬ 


ential  distribution  is  also  a  limiting  case  of  Ila. 
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The  properties  exhibited  by  this  very  general  class  of  functions 
are  made  use  of  in  the  simpler  of  the  two  Bayesian-dynamic  models  de¬ 
veloped  subsequently.  In  the  example  presented,  the  simplest  a  priori 
assumption  is  made,  namely  that  the  probability  of  survival  to  time  t 
is  exponential.  However,  most  models  could  just  as  well  involve  more 
general  functions  which  are  members  of  this  class. 

In  subsequent  chapters  we  develop  these  models. 


CHAPTER  VI:  INVENTORY  MODELS  WITH  MARKOV  DEMAND 


A.  Introduction 

In  this  chapter,  we  present  some  inventory  models  which  are 
useful  for  decision-making  regarding  ordering  and  disposal  activities 
involving  a  single  item.  Policies  considered  are  of  the  well-known 
(s,S)  type.  For  ordering,  whenever  the  stock  level  falls  below  a 
certain  level  s,  enough  new  stock  is  ordered  to  bring  the  level  up 
to  another  prespecified  level  S,  if  the  stock  level  exceeds  s, 
no  order  is  placed.  As  to  disposal,  the  policy  operates  as  follows: 

If  the  stock  level  exceeds  a  certain  level  D,  a  disposal  takes 
place  in  order  to  bring  the  stock  level  down  to  a  new  level  d;  when 
the  level  is  short  of  D,  no  disposal  is  made. 

In  Section  2,  a  demand  process  is  formulated.  We  assume  that 
the  system  which  generates  demand  can  be  in  any  of  several  states  in 
a  given  time  interval.  Each  state  has  its  own  demand  pattern.  Hence, 
the  demand  in  each  period,  which  is  a  random  variable,  may  or  may  not 
be  identically  distributed  in  successive  periods.  At  the  end  of  each 
time  interval,  we  have  data  on  the  demand  for  that  interval^  based  on 
these  data  we  then  efficiently  forecast  the  demand  for  the  future  by 
means  of  Bayes*  formula. 

In  Section  3,  a  general  inventory  model  which  is  embedded  in  the 
demand  process  described  in  the  previous  section  is  presented  with  a 
computation  procedure.  A  specific  example  of  such  a  model  is  given 
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witfa  some  calculated  results.  Finally,  a  very  special  case  of  the 
general  model  is  discussed.  This  last  model  is  distinguished  by  the 
property  that  when  the  system  has  made  a  transition  to  its  terminal 
state,  which  is  characterized  by  zero  demand,  it  stays  there.  In 
Section  U ,  two  schemes  which  can  be  used  for  improving  the  accuracy 
of  predicting  the  demand  for  the  future  are  introduced. 

B.  The  Demand  Process 

1  .  Markov  Demand  Generation 

Of  course,  the  system  which  underlies  the  demand  generation 
of  even  a  single  item  in  the  vast  Navy  supply  system  is  quite  complex. 
Fortunately,  a  Markov  process  provides  a  manageable  and  quite  general 
mathematical  model  for  our  analytical  study. 

The  Kasic  concepts  of  the  Markov  processes  are  those  of  "states" 
of  a  system  and  state  "transition."  We  say  a  system  is  in  a  certain 
state  when  it  is  completely  describable  by  the  values  of  parameters 
which  define  that  state.  A  system  is  said  to  have  made  a  transition 
from  one  state  to  another  when  the  parameters  which  describe  it  have 
changed  from  the  values  specified  for  one  state  to  those  for  another. 

Consider  a  demand  generating  system  which  can  be  in  one  of  a 
finite  number  of  states,  where  r  =  l,2,...,n,  at  any  time. 

Each  state,  in  general,  is  assumed  to  have  a  different  demand  pattern 
which  is  characterized  by  a  parameter  such  as  mean  demand  rate.  This 
system  makes  state  transitions  according  to  a  certain  transition 
probability  matrix  j|P  ||.  Its  typical  element  stands  for 


|cQ| 
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the  probability  of  the  system  making  a  transition  from  S  to  S  . 

IT  S 

Furthermore,  we  use  a  vector  («  )  to  denote  the  a  priori  state 
probability  distribution.  Its  typical  element  jc  is  the  proba¬ 
bility  that  the  system  is  initially  in  S  and  Z  jt  =1.  Let  g  (x) 

it  r  r  it 

denote  the  conditional  distribution  of  demand  X  given  that  the  sys¬ 
tem  is  in  S  . 

r 

The  above  finite  state  discrete  time  process  is  the  scheme  in¬ 
corporated  in  the  subsequent  inventory  models  to  represent  the  de¬ 
mand  generation  for  a  single  item.  This  process  can  also  readily 
handle  the  situation  where  demands  for  several  items  are  linked 
through  underlying  states.  All  one  must  do  is  specify  a  condi¬ 
tional  joint  distribution  of  item  demands  for  each  state.  One  of 
the  models  to  be  discussed  later  makes  use  of  this  notion  of  linked 
demands . 

2.  Prediction  of  Demand  by  Bayes'  Formula 

Suppose  we  have  observations  on  demand  X  in  this  period 
and  we  wish  to  estimate  the  demand  for  the  next  period.  How  shall 
we  proceed?  One  method  of  prediction  is  to  calculate  a  posteriori 
state  probabilities  by  means  of  Bayes  formula.  Then  combining  these 
state  probabilities  with  the  transition  probabilities  provides  a  basis 
for  estimation. 

We  are  given  the  a  priori  star e  probabilities  «r=Pr  [s=sr]  f°r 
r  =  1,2, ...,n  and  we  wish  to  calculate  the  a  posteriori  state 


0*1 


□  « 
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probabilities  =  Pr  £s  =  Sr  j  X  =  xj  ,  for  r  =  1,2, ... ,n,  after 
observing  X  =  x.  Using  the  definition  of  conditional  probability 

=  Pr  [s  =  Sr  |  X  =  x]  , 

li: 


(6.1) 


Pr  S  =  S  and  X 
1  r 


-*] 


Pr  £x  =  xj 

The  numerator  of  (6.1)  can  be  evaluated  as  follows: 

(6.2)  Pr|s  =  Sr  and  X  =  xj  =  Pr^X  =  x|S  =  Srj.Pr|s  =  Sr]  , 

=  gr  W  •  *r. 

To  evaluate  the  denominator  we  note  that  the  possible  outcomes  leading 
to  X  =  x  are  (S^x),  (S2,x),  .  .  .  (Sn,x)  Hence 

(6.3)  Pr^X  =  x]  =  ^  Pr[x  =  x  and  S^j  , 


=  2  g  (x)  .  jt  . 

j=l  J  J 


From  (6.2)  and  (6.3),  we  can  then c alculate  (6.1). 
(6.4) 


"r  n 


gr  (x)  -  *r 


gj  (x)  • 

Making  use  of  knowledge  about  the  a  posteriori  state  probability 
and  the  state  transition  probability  matrix  ||Prs||,  we  then  obtain 
the  a_  priori  state  probability  (n^')  for  the  next  period. 
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(6.5) 


n 

it"-  2  P 
s  r=i  rs 


« ", 

r' 


A  •  8r  <»>  '  »r 
j5i  h  <x)  •  -j 


(Kg')  gives  us  a  basis  for  predicting  the  future  states  of  the  system 
on  the  basis  of  all  observable  demands  up  to  the  present,  by  recursive 
application  of  (6.5)  after  each  time  period.  Furthermore,  it  also 
makes  use  of  available  a  priori  information  regarding  the  future  states 
of  the  system.  For  instance,  so  called  program  data  such  as  the 
expected  numbers  of  flight  hours  and  the  expected  obsolescence  rate 
may  be  translated  into  the  form  of  the  a  priori  state  probabilities 
(«r)  and  the  state  transition  probabilities  ||prsl|*  Pn  this  way, 
the  empirical  evidence  can  be  supplemented  by  the  judgment  of  the 
operating  people  to  obtain  the  most  efficient  prediction  results. 

C .  General  Markov  Inventory  Model 
1.  Formulation  of  the  Model 

The  problem  considered  here  is  that  of  stocking  a  supply  of 
one  item  to  meet  an  uncertain  demand  which  is  assumed  to  be  generated 
by  a  Harkov  process  described  in  the  previous  section.  Various  costs 
associated  with  oversupply  and  under supply  are  assumed  to  be  operative. 
Orders  are  placed  at  the  beginning  of  each  time  period  of  equal  length. 
The  orders  are  assumed  to  be  fulfilled  either  immediately  or  are 


|cQ 
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delivered  one  time  interval  later.  After  the  order  has  been  placed,  a 
demand  is  made.  This  demand  is  satisfied  from  the  existing  inventory, 
with  excess  demand  leading  to  a  penalty  cost.  Unfilled  demand  can  be 
assumed  to  be  either  lost  or  backlogged*. 

Let  us  first  introduce  the  following  notation: 

\ji  =  policy 

a  =  action  according  to  policy  presumably  a  =  t(Y,  (jt)) 

A  =  the  set  of  admissible  actions,  aeA 

Y  =  inventory  vector,  describing  the  inventory  status  at  the 
beginning  of  a  period 
X  =  demand  vector 
Jt  =  a  priori  state  probability 

L(Y,x,a)  =  expected  current  cost  function  with  initial  inventory 
Y,  demand  x  and  taking  action  a. 

C^(Y,  (rc^_))  =  expected  total  cost  of  choosing  a  according  to 
and  process  beginning  with  Y  and  (jt^)  . 

The  optimization  problem  involves  finding  to  minimize 
Cy  (Y,  (itp)  for  all  (Y,  (it^)).  We  note  the  following  recursive 
relation: 

(6.6)  Ct(Y,(*r))  =  gZx  jtsgs(x)jL(Y,x,a)  +  P  Cy  (Y‘,  (*”))]  , 

^As  long  as  this  assumption  is  maintained,  the  computation  of  opti¬ 
mal  policies  for  any  constant  lead  time  problem  is  still  manageable. 

As  to  the  inventory  problem  with  a  stochastic  lead  time,  the  literature 
is  rather  scarce  on  this  subject  except  under  some  simplified  assumptions 
clearly,  further  exploration  is  needed  in  this  area  of  inventory  study. 


Y1  =  a  vector  describing  terminating  inventory  status  for 
this  period  (or  initial  inventory  status  for  next 
period) a  It  assumes  the  following  functional  rela¬ 
tione  Y'  =  h  (Y,  x,  a) 

jt"  =  a  priori  state  probability  for  next  period.  This  was 

derived  in  the  previous  section,  as  a  function  of 

(jt  )  and  x . 
s' 

p  =  a  discount  factor,  0  <  p  <  1.  The  introduction  of 
such  a  discount  factor  prevents  infinite  costs  from 
entering. 

Following  the  well  known  optimality  principle,  the  optimum  total 
cost  function  may  be  stated  as  follows; 

(6.7)  Cv(Y,0tr))  =  Min  Z^  «sgg(x)  [l(Y,x,«)  +  p  C^(Y',  (jtJJ))] 

If  t  is  an  optimum  policy. 

This  functional  equation  can  be  solved  by  the  following  iterative 
method;  Select  an  arbitrary  set  of  starting  policies  and  cost  func¬ 
tions,  say  H,o  and  Kq  (Y,  (nf)).  Then  recursively  calculate  sequen¬ 
ces  of  policies  and  cost  functions  ^Kn(Y, (nr) )^  such  that 

(6.3)  Vi(Y,(*r))=Min  sEx  *8gs(x)  [l(Y,x,o)  +  P  Kn(h(Y,x,C«),  («”))] 

and  i|r  (Y,  (it  ))  is  the  minimizing  a  for  each  (Y,  (jt  )).  When 
ht  i  r  r 

these  sequences  converge,  the  limits  of  the  sequences  are  the  solution 
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to  (6o7)o  In  our  numerical  study  of  this  model,  it  was  found  that 
the  convergence  of  the  cost  function  sequence  could  be  pain¬ 

fully  slow,  expecially  when  the  values  of  starting  cost  functions 
are  poor  approximations  of  the  true  cost  functions  or  the  discount 
factor  ,p  is  close  to  unity.  However,  it  was  found  that  the 
policy  sequence  converges  much  more  rapidly. 

2.  Computational  Short-cut 


We  note  that  in  any  minimization  or  maximization  problem 
of  this  type  what  matters  is  the  pattern  of  differences  among  the 
relevant  variables,  rather  than  absolute  levels.  Making  use  of  this 
idea,  we  are  then  able  to  reduce  the  number  of  iterations  necessary 
for  the  sequence  ^  K  to  converge. 

Let  (Y,(*r))  be  any  particular  value  of  (Y,  (it^)),  and  define 

(6  =  9)  C+  (Y,(nr))  =  <Y,(*r»  -  Cj  (?,(*r», 

where  i|r  is  optimum  policy  and  corresponding  minimum  cost. 

Substituting  (6.9)  into  (6.7),  we  have 

(6.10)  Ct(Y,(«r))  =  M^n  sLx  *sgg(x)  [l(Y,x,0c)  +  P  (h(Y,x,o)  ,  («”) )] 


-  (1-P)  (Y,(«r)). 


This  subtraction  of  a  constant  will  in  no  way  affect  the  mini- 

°  o 

mization.  Hence,  we  set  C^(Y,  («r))  =  0.  We  then  have  another  expect¬ 
ed  total  cost  function,  involving  only  the  cost  differences. 
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(6 •  11)  (Y,(itr))  =  M^n  s^x«sSs(x)  [l(Y,x,Cx)  +  P  (h(Y,x,a)  ,  (*£))] 

The  new  iterative  procedure  is  exactly  the  same  as  before  except 
after  each  iteration  the  following  subtraction  is  made  to  set 

Vi&<*r»  "  ° 

(612)  K^i (Y,  (if r} )  =  Min  Z*  it^s )  [l(Y,x,cO  +  P  Kn(h(Y,x*,a) ,  (*;;))] 

'MinSjX  7ts8s(x)  |_L(Y,x,a)  +PKn(h(Y,x)a),(np)] 

When  convergence  of  the  sequence  is  satisfactory,  the  last  sub- 

o  o 

tractive  constant  is  divided  by  (1-p)  to  obtain  C^(Y,  . 

o  o  Min  S  k  g  (x)  fL(?,x,a)  +  P  K  (h(Y,x,a),  (it"))] 

(6  13)  C  (Y,  (it  )  )  =  -3-^-x-S  5 - ^ - S - 

V  1-P 

2 

This  is  a  quick  way  of  summing  the  infinite  series  (1  +  P+  p  +. .  .) 
Finally,  (Y,  («r) )  for  other  (Y,  (Jtp)  are  calculated  by  means  of 
(6,9) 

By  adopting  this  short-cut  method,  a  substantial  reduction  in  com¬ 
putation  time  was  observed. 

3.  Numerical  Examples 

In  this  numerical  example,  a  demand- gener at ing  system  consist¬ 
ing  of  two  states,  one  having  a  high  demand  rate,  and  the  other  a  low 
demand  rate,  is  considered.  Orders  are  made  at  the  beginning  of  each 
regularly  spaced  period  and  a  delivery  lag  of  one  period  is  considered. 
Any  unfilled  demand  for  current  period  becomes  demand  for  next  period. 

The  conditional  distribution  of  demand  is  assumed  to  be  Poisson. 


Definition  of  symbols  used  and  parameter  values  assumed  are  given 
below: 

y  =  initial  inventory  level.  Both  positive  and  negative  values 
are  allowed.  If  y  assumes  a  negative  value,  it  means  that 
there  are  unfilled  demands. 

a  =  desired  inventory  level.  The  difference  a-  y  represents 
the  optimal  order  quantity,  if  it  is  positive.  If  the  dif¬ 
ference  is  negative,  it  represents  the  optimal  disposal  quan¬ 
tity.  In  the  first  example,  we  are  considering  a  >  max  (y,0), 
i.e.,  no  disposal  activities  are  allowed.  Of  course,  it  pre¬ 
sents  no  problem  if  one  wishes  to  allow  disposal  activities 
as  in  the  second  example)  all  that  needs  to  be  done  is  to  set 
the  domain  of  a  to  be  a  >  0.  If  we  allow  a  to  assume 
negative  values  it  means  that  the  stock  can  be  returned  to  the 
suppliers  and  credit  received  for  it. 

x  =  demand  per  period 

n  =  a  priori  probability  that  the  system  is  in  state  1.  Then 

(1-jt)  is  a  priori  probability  that  the  system  is  in  state  2. 
g^(x)  =  probability  distribution  of  demand  when  the  system  is  in  state  1 

-2  2X 
6  71 

g£(x)  =  probability  distribution  of  demand  when  the  system  is  in  state  2 

-0.4  (0.4)x 
e  _ 1 
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probability  distribution  of  x  weighted  by  a  priori  state 
probabilities 

hj (x)  =  x  gx(x)  +  (1-x)  g2(x) 
unit  storage  cost,  0.5 
unit  out-of-stock  cost,  5.0 
fixed  order  cost,  1.0 
unit  order  cost,  0.5 

transition  probability  matrix,  r,s,  =1,2 


next  period. 


Pn  *  Mx>  +  p2i  t1-*)  gz(x) 

^  ~  «  g1(x)  +  (1  -  n)  g2(x) 

p  =  discount  factor,  0.99 

With  the  above  notation,  we  now  can  proceed  to  specify  the  expected 
total  cost  function  which  is  composed  of  the  expected  current  costs  and 
discounted  future  cost. 

Expected  holding  and  shortage  costs,  L^(y,jt),  to  be  charged  dur¬ 
ing  the  current  period,  assuming  that  an  order  will  not  be  delivered 
until  next  period  is: 
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(6.14)  Ll(y,n)  =« 


'dl  x=0(y'X)  h«  (X)  +  d 


<*2  (-y)  +  d  f, 

*■  1  x=0 


Suppose  the  ordering  cost,  L2(y,  a), 
placed,  it  can  be  expressed  as  follows: 


L2(y 


j  a)  ={ 


d3  +  d4  (a  -  y) 


,  ^  (x-y)  h  (x)  y  >  0 

■x=y+l  n 

x  (x)  y  <  0 

is  charged  when  orders  are 

a  >  y 

a  =  y 


Then  the  expected  current  cost  function  is: 

Hy,  *>  a)  =  L1(y,  Jt)  +  L2(y,  u) 

Let  (y,>:)  represent  the  total  expected  cost  if  the  provisioning  is 

done  optimally. 


(6.15)  C.  (y,Jt)=  Min  |L(y,it,a)  +  P  JL  h 
V  a>max(y,0)  x  u  n 


00 


(a-y,Jt" 


,1 


The  solution  to  (6.15)  was  obtained  by  means  of  the  iterative  pro¬ 
cedure  described  in  the  previous  section  and  is  presented  in  Table  1. 
Note  that  the  condition  a>max(y,0)  restricts  the  policy  from  carrying 
out  disposal  activities. 

Another  set  of  calculations  was  performed  by  introducing  dis¬ 
posal  activities  in  the  above  problem  and  by  reducing  the  unit  stock¬ 
out  cost.  (Rest  of  the  parameters  were  unchanged.)  More  specifically, 
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TABLE  1:  OPTIMUM  POLICY  TABLE  WITHOUT  DISPOSAL 


y  \  o.i 


a  priori  state  probability 

0.3  0.4  0.5  0.6  0.7  0.8 


/ 

6 

/ 

6 

/ 

6 

/ 

6 

/ 

6 

7 

6 

7 

6 

7 

6 

7 

6 

5 

5 

5 

5 

5 

5 

5 

5 

5 

4 

4 

4 

4 

4 

4 

6 

6 

6 

3 

5 

5 

5 

5 

5 

6 

6 

6 

4 

5 

5 

5 

5 

5 

6 

6 

6 

4 

5 

5 

5 

5 

5 

6 

6 

6 

4 

5 

5 

5 

5 

C 

J 

6 

6 

6 

4 

5 

5 

5 

5 

5 

6 

6 

6 

4 

5 

5 

5 

5 

5 

6 

6 

6 

For  a  given  a  priori  state  probability  and  initial  inventory  level,  we 
read  off  an  appropriate  entry  in  the  above  table.  The  difference  between 
this  entry  and  the  initial  inventory  level  is  the  optimum  order  quantity. 
When  the  difference  is  zero,  no  action  should  be  taken.  For  all  y  <  -2, 
the  optimum  policy  repeats  itself. 
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we  set  the  domain  of  minimization  to  a  >  0  in  (6.15)  and  set  =  25. 
Results  are  tabulated  in  Table  2. 


TABLE  2:  OPTIMUM  POLICY  TABLE  WITH  DISPOSAL 


a  priori  state  probability 
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has  occurred  when  the  system  has  entered  the  zero  demand 
state. 

2.  It  illustrates  how  those  failure  or  mortality  distributions 
with  the  property  of  "conditional  reproducibility"  which 
have  been  discussed  in  an  earlier  chapter,  enter  into  the 
general  Markov  inventory  model.  With  the  introduction  of 
such  a  distribution  into  the  model,  it  is  appropriate  to 
look  for  optimum  policies  which  depend  on  the  inventory 
level  and  on  some  function  of  the  past  observable  demands 
but  not  explicitly  dependent  on  time. 

Consider  a  system  with  two  underlying  states.  When  the  system  is  in 
State  1,  the  item  demand  follows  the  conditional  distribution  g^(x) 
for  x=0,  1,2,...  When  it  is  in  State  2,  the  conditional  demand 
function  g2(x)  is  defined  to  be 


g200  = 


1  for  x  =  0 

.0  for  >  0 

Then,  in  the  language  of  the  previous  section,  the  probability 
distribution  of  X  weighted  by  a  priori  state  probabilities  is 

for  x  =  0 

.«  gx(x)  x  >  0 


h  (x)  = 


*  gx(x)  +  (l-«) 


Suppose  transitions  are  made  according  to  the  following  transition 
probabilit ies 
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Pr(t  <  T  <  t+dt|T  >  t)  =  The  probability  that  the  system  will  terminate 

in  an  additional  time  interval  dt  given  that  it  did  not 

terminate  before  t. 

-kt 

Suppose  we  assume  f(t)  =  ke  .  Then 


Pr (T  <  t)  =  1  -  e 


Pr(t  <  T  <  t  +  dt)T  >  t)  *  1  -  e'k(dt) 

Now  the  probability  that  the  system  will  survive  the  additional 
time  interval  dt  given  that  it  has  survived  up  to  t  is 

1  -  Pr(t  <  T  <  t  +  dtjT  >  t)  =  e"k^dt^ 

If  we  set  dt  to  be  unit  time,  we  have  the  desired  probability 


D.  Some  Extensions  of  the  General  Model 
1.  A  Model  with  Linked  Demands 

The  general  model  is  extended  to  handle  linked  demands,  in 
which  linkage  is  based  on  underlying  state  S^.  for  r  =  l,2,„..,n, 
in  the  following  manner.  With  each  state  we  associate  a  joint  density 
function  of  demands  such  that  when  the  system  is  in,  say,  Sg  the  re¬ 
sulting  dmeands  are  considered  to  be  observations  based  on 

8s  (xl,x2  ’  °  *  »xx  >  ' 

If  the  information  such  as  that  mentioned  above  can  be  incorporated 


into  our  model,  it  will  contribute  to  the  statistical  determination  of 
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where  the  system  is  in  any  one  period.  For  instance,  if  we  know 
_a  priori  that  demands  for  items  A  and  B  occur  together  in  state 
i  but  in  state  j  there  is  demand  only  for  item  A,  and  then  if 
there  are  observable  demands  for  item  B,  we  will  conclude  that 
the  system  is  most  likely  to  be  in  state  i. 

With  the  introduction  of  linked  demands,  the  calculations  of 
the  a  posteriori  state  probability  (jt^)  and  the  new  a  priori  state 
probability  (jt”)  (as  in  (6A)  and  (6.5)  will  be  modified  as  follows: 

8r(W  '  *  Mxi)Jtr 

*  '  =  -  > 

r  n 

jfl 

n 

r=l  Prs  8r(V-V"r 
**  =  _  . 

jfx  8j 

As  to  optimization,  it  will  be  carried  out  for  one  item  at  a 
time.  Necessary  modifications  are  quite  straightforward. 

Suppose  we  are  interested  in  optimization  of  the  ith  item. 

The  expected  total  cost  function  can  be  stated  as  follows:  (Note  all 
the  variables  have  the  same  meaning  as  before  except  now  they  refer 
to  the  ith  item) . 

(6.16)  Vyi,{nr}  }  =  a\n  5tsms(xi)fL(yi'xi’ai)  +PVyMns}) 
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where  mg(x^)  is  the  marginal  distribution  of  x^ 

00  00 

ms (xi>  =  /•••••  f  Ss(xr--xp  d  xr--dxi-i  dxi+i---dxi- 

O  o 

Again,  this  functional  equation  can  be  solved  by  the  same  iterative 
procedure  described  earlier. 

2 .  A  Model  with  Non-demand  Ancillary  Information 

The  general  model  can  also  include  non-demand  ancillary 
information  which  is  indicative  of  state  of  the  system.  This  will  be 
illustrated  with  the  following  example. 

Consider  a  three-state  system  in  which  the  first  two  states  have 
the  same  demand  pattern  but  different  state  transition  probabilities. 

The  third  state  follows  its  own  demand  pattern  and  transition  proba¬ 
bilities  different  from  the  first  two  states. 

In  such  a  case,  since  observable  demands  alone  cannot  distinguish 
the  first  state  from  the  second  state,  observations  on  non-demand 
variables  which  have  probability  distributions  dependent  on  the  state 
will  be  useful  for  sharper  statistical  discrimination. 

Consider  the  following  example:  the  ancillary  information  avail¬ 
able  is  about  conditional  distribution  of  a  non-demand  variable  T. 

Pr  1  T  =  tjS  =  S  |  ,  where  t=l,2,3 

•  r  and  r=  1 ,  2 ,  3 . 

This  distribution  may  be  assumed  to  look  like  the  following  table: 
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u 
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Each  row  defines  a  distribution  of  T  in  the  corresponding  state. 

The  variable  T  can  be  interpreted  as  being  a  more  or  less  accurate 

indicator  of  the  demand  state.  Such  information  is  believed  to 

be  readily  available  within  the  Navy  Supply  System  at  this  time. 

Next,  given  the  conditional  probability  Prj T  =  1 1 S  =  j,  it 

is  of  interest  to  know  what  is  the  a  posteriori  probability 

Pr  S  =  S  It  =  t  i  which  is  the  conditional  distribution  of  the  states 
r 1 

of  the  system  given  the  ancillary  information  T  =  t.  This  can  be 
calculated  readily  from  the  definition  of  conditional  probability. 


Prob.  of  T 


\  1  1 
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3 

S1 
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Pr;s  =  s  It 
l  r 1 


Let  Pr|s  =  Sr|T 


Pr'S  =  S  and  T  =  t 
1 _ r _ 

Pr|  T  =  tj 


Pr  T  =  t| 

S  =  s  1  . 
r 

Pr^S 

«  S  1 
r 

£  Pr |T  = 
i 

=  t  and  S 

=  s.n 
1 

J 

Pr  T  =  t| 

S  =  S  |  . 
r  I 

Pr  S 

"  Sr  ] 

Z  Pr  T  = 
i 

t|s  =  S' 
J 

.  Pr 

'S  =  Sjj 

t]  =  Kr  (T) 
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«r(T)  will  then  replace  the  a  priori  state  probabilities  n  (p.58) 

in  the  calculations  of  the  a  posteriori  state  probabilities 
f  1 

Pr j S  =  Sr|x  =  x  j y  and  the  a_ priori  state  probabilities  of  next  period 

p  Pr  S  =  S  lx  =  xj. 
rs  r1  J 

Optimization  follows  exactly  the  same  steps  as  in  the  general 
model . 
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CHAPTER  VII 


SOME  RESULTS  OF  A  STATISTICAL  STUDY  OF  OBSOLESCENCE  DATA 

A<  Introduct Ion 

In  order  to  generate  optimal  inventory  decision  rules  from  some  of 
the  models  described  earlier  in  this  report  (see  Chapters  IV  and  VI), 
it  is  necessary  to  assume  some  a  priori  probability  distribution  of 
the  occurrence  of  obsolescence  of  an  item.  While  it  may  be  true  that 
incorrect  assumptions  regarding  the  probability  distribution  of  obsoles 
cence  is  self-correct ing  in  the  long  run,  as  we  accumulate  more  experi¬ 
ence  with  the  item,  the  fact  that  we  are  dealing  with  a  process  which 
is  only  finite  in  duration  suggests  that  the  more  we  know  a  priori 
the  more  likely  are  the  decisions  to  be  optimal.  We  are  also  inter¬ 
ested  in  enhancing  our  understanding  of  the  factors  which  are  associa¬ 
ted  with  the  process  of  obsolescence,  i.e.,  we  wish  to  determine 
causal  factors  in  the  process  which  set  an  obsolescence  date;  if  this 
proves  impossible,  we  wish  to  be  able  to  at  least  isolate  those  fact¬ 
ors  which  are  related  to  the  process.  Therefore,  data  bearing  on 
life-spans  of  a  sample  of  items  (mainly  airframe  and  engine  parts), 
along  with  a  variety  of  economic  and  technical  data  about  each  of  the 
items  was  examimed. 

The  items  were  first  sorted  according  to  Federal  Stock  Code  (FSC) 
and  Technical  Supply  Maintenance  Code  (TSMC).  Within  these  groups  we 
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studied:  (1)  the  frequency  distribution  of  time-to-obsolescence,  and 
(2)  the  degree  of  associations  between  some  economic  and  technical 
variables  and  life-spans. 

We  learned  that  a  certain  type  of  the  failure  distributions  appears 
to  describe  some  of  the  observed  data  adequately  and  that  there  is  a 
statistically  significant  relation  between  the  life  spans  of  the  items 
in  some  groups  and  certain  factors  associated  with  the  items.  However, 
the  data  examined  indicated  that  there  was  no  single  frequency  distri¬ 
bution  that  may  be  usefully  employed  across  the  board  to  describe  the 
pattern  of  the  occurrence  of  obsolescence.  This  emphasizes  the 
difficulty  of  predicting  when  obsolescence  will  take  place. 

B.  Data  Analyzed 

Before  entering  into  the  discussion  of  the  results  of  the  statisti¬ 
cal  analysis,  we  shall  describe  in  this  section  the  data  used  in  this 
study  and  their  limitations. 

The  data  analyzed  were  obtained  primarily  from  the  Aviation  Supply 
Office  (ASO).  Initially  a  list  cf  about  8,000  items  which  were  declared 
obsolete  during  the  first  half  of  1961  was  prepared  by  ASO  for  this 
study.  Unfortunately,  this  list  was  not  very  useful  because  most  of  the 
information  relating  to  these  items  was  no  longer  available  when  the 
data  collection  started  in  September  1961.  An  alternative  list  of  more 
recent  obsolescence  items  was  then  compiled  from  the  Stock  Number  Data 
Section  (SNDS)  catalogs  issued  in  August  and  September  of  1961.  These 


items  were  declared  obsolete  during  the  13-week  period  ending  with  the  dat 
of  SNDS. 


The  additional  sources  of  information  pertaining  to  the  behavior  of 
the  obsolescence  items  in  SNDS  were  found  in  the  following  documents:  The 
Retention  and  Disposal  Listing,  the  Parts  List  Catalog,  the  Purchase 
Order  History  Cards  (also  called  7-30  cards  or  On-Order  Cards),  the 
Federal  Supply  Classification  Numeric  Index  of  Classes,  the  Naval  Air¬ 
craft  Maintenance  Program  Glossary,  and  some  specially  prepared  lists 
matching  Navy  contract  numbers  to  their  dates  of  issues. 

The  data  obtained  from  the  above  documents  are  by  no  means  complete. 
For  instance:  the  Retention  and  Disposal  Listing  includes  only  data  on 
items  with  surplus  stock  prior  to  the  declaration  of  obsolescence;  the 
parts  list  catalogs  do  not  list  all  the  applications  for  some  of  the 
items;  between  25 %  and  75%  of  the  On-Order  cards  for  obsolete  items  were 
missing  from  the  files.  The  variation  in  missing  cards  seems,  in  large 
part,  to  reflect  the  efficiency  and/or  prudence  of  the  clerks  responsible 
for  them.  There  is  a  policy  of  purging  cards  from  the  files  of  those 
items  which  have  not  had  purchase  orders  in  the  preceding  two-year  period 
Hence,  some  relevant  information  for  those  items  which  have  not  had  a 
regularly  recurring  pattern  of  buys  is  missing. 

Despite  the  incomplete  nature  of  most  of  the  data,  there  was  On- 
Order  card  data  for  about  4,000  items.  For  all  these  4,000  items  there 
is  supplementary  information  from  the  Part  List  and  the  SNDS.  For  ap¬ 
proximately  1,500  of  these  items  there  is  also  supplementary  information 
from  the  Retention  and  Disposal  Listing  for  1961.  In  addition,  we  were 
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able  to  obtain  the  relevant  information  for  an  additional  1,000  obsolete 
items  from  the  Douglas  Aircraft  Company  at  El  Segundo,  California. 

Information  Contained  in  the  Data 

The  following  list  contains  the  descriptive  and  quantitative  vari¬ 
ables  available  for  the  study,  their  limitations,  and  their  source 
documents.  The  variables  are  randomly  ordered  with  no  particular 
thought  to  the  relative  importance  of  any  one  of  them. 

Federal  Stock  Number:  (FSN)  A  17  digit  alphanumeric  interservice 
part  number  by  which  all  items  are  classified  in  all  the  catalogs 
used  for  the  study.  The  third  through  the  sixth  digits  are  numerical 
and  denote  the  Federal  Supply  Class.  The  fourteenth  through  the 
seventeenth  digits  are  alphabetic  and  denote  the  Technical  Supply 
Maintenance  Code,  the  last  three  digits  of  which  represent  the 
manufacturer.  The  other  digits  are  not  significant  for  this  study. 

A  Federal  Supply  Class  incorporates  those  items  which  support 
similar  functions,  although  the  particular  items  are  dissimilar, 
e.g.,  4120  (Fire  Fighting  Equipment)  includes  axes,  ladders,  and 
firehouse  carts.  Some  classes,  however,  are  restricted  to  similar 
parts  only,  e.g.,  6240  (Light  Bulbs  and  Lamps). 

Unit  Price:  The  average  price  paid  for  a  given  item  in  a  given 
contract.  Both  setup  arid  manufacturing  costs  are  imputed  in  the 
unit  price  for  each  separate  contract.  Hence  the  unit  price  varies 
according  to  changing  labor  and  setup  costs  associated  with  the 
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different  quantities  purchased  at  different  points  in  time. 

The  unit  price  is  sometimes  listed  in  the  SNDS  and  the  On- 
Order  Cards;  for  the  On-Order  Card,  several  different  prices 
may  be  noted.  A  single  price  always  appears  in  the  Retention 
and  Disposal  Listing  and  the  Parts  List  Catalog.  The  choice 
of  price  for  the  Parts  List  is  not  clear. 

Nomenclature:  A  descriptive  name  for  a  part,  which  is  designed 
to  group  an  item  by  its  physical  characteristics  and/or  function. 
For  the  study  the  groups  may  be  overly  definitive,  i.e.,  shims, 
washers,  and  spacers  may  be  sufficiently  similar  as  to  con¬ 
stitute  a  single  group.  All  the  catalogs  and  cards  indicate  the 
nomenclature  of  an  item. 

Source  Code:  A  code  denoting  the  origin  of  a  part,  e.g., 
internal  manufacturer,  inter service  transfer,  or  commercial 
contractor  is  found  both  in  the  SNDS  and  the  Parts  List. 

BUSANDA  Change  Code:  A  code  found  in  the  SNDS  indicating  the 
reason  that  an  item  was  declared  obsolete  if  preceded  by  a  "Q". 

Life-Span:  The  number  of  months  between  initial  BUWEPS  procure¬ 
ment  of  an  item  and  the  declaration  of  obsolescence.  This  number 
may  be  computed  from  the  data  of  the  BUWEPS  contract  (prefixed 
"NOAS")  listed  first  on  the  On-Order  Card.  If  the  date  is  absent 
from  the  card,  it  may  be  obtained  from  a  list  of  BUWEPS  contracts. 
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Since  the  early  On-Order  Cards  for  many  Items  are  often  missing, 
two  problems  arise.  First,  no  NOAS  contract  number  may  appear 
on  the  card,  indicating  that  the  first  date  on  the  card  is  for 
a  re-order  subsequent  to  initial  procurement.  Second,  an  item 
with  several  applications  will  have  had  several  initial  procure¬ 
ments,  i.e.,  NOAS-coded  contracts,  one  for  each  application. 
Hence  the  first  NOAS-coded  contract  noted  on  a  card  may  not  be 
the  first  for  that  item.  With  either  problem,  the  life-span  for 
an  item  will  be  understated.  Items  with  the  first  problem  may 
be  flagged  to  indicate  that  the  life-span  is  understated.  Those 
with  the  second  problem  cannot  be  so  differentiated. 

Quantity  on  Hand:  That  quantity  of  each  item  remaining  at  all 
Navy  Supply  Depots  at  the  declaration  of  obsolescence.  This 
figure  is  not  available  from  any  source.  One  approximation  for 
each  item  is  the  quantity  on  hand  at  the  last  annual  running  of 
the  Retention  and  Disposal  program.  The  program  computes  an 
optimal  maximum  retention  quantity  for  each  item.  Any  quantity 
in  stock  greater  than  this  maximum  is  termed  surplus.  Provided 
that  a  surplus  existed  at  the  time  of  the  running,  the  item  will 
appear  on  the  Retention  and  Disposal  listing,  giving  the  maximum 
retention  quantity  and  the  quantity  on  hand.  The  most  recent 
running  of  the  program  took  place  about  six  months  prior  to  the 
publications  of  the  SNDS  catalogs  used  for  this  study.  Hence, 
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inclusion  in  the  study  only  of  items  with  "on- hand"  date  would  give 
a  bias  towards  items  prone  to  surplus  and  discount  all  quantities  of 
items  withdrawn  subsequent  to  the  program,  but  prior  to  the  declara¬ 
tion  of  obsolescence.  Inclusion  of  the  balance  of  the  items  re¬ 
quires  that  a  zero  on-hand  quantity  must  be  assigned  to  those 
items,  which  modifies  the  bias  to  those  items  which  either  had  a 
surplus  or  whose  stocks  were  exhausted  at  the  precise  moment  that 
the  declaration  of  obsolescence  took  place.  It  is  not  possible  to 
isolate  those  items  whose  stocks  had  been  exhausted  prior  to  the 
declaration,  and  were  unable  to  satisfy  issue  requests. 

Number  of  Purchases:  The  number  of  times  that  re-orders  were  made 

is  noted  on  the  On-Order  Cards- 

Total  Purchases:  The  sum  of  all  the  purchases  made  is  found  by 
adding  the  entries  on  the  On-Order  Cards. 

Maintenance  Per  Cent:  The  percentage  of  anticipated  replacements 
during  maintenance  is  available  from  the  Parts  List. 

Overhaul  Per  Cent:  The  percentage  of  anticipated  replacements 
during  overhauls  is  also  given  on  the  Parts  List. 

Application  Codes:  A  code  giving  either  the  end  item,  i.e.,  the  plane, 
on  which  the  item  goes,  or  the  assembly  of  which  the  item  is  a  replace¬ 
able  part.  This  code,  found  in  the  Parts  List,  gives  a  measure  of  the 


interchangeability  of  a  part  within  one  plane,  or  among  many  planes. 
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Unlts  Pet  Application:  The  quantity  of  a  particular  item  needed 
to  perform  the  desired  function  is  found  in  the  Parts  List.  It 
is  possible,  although  difficult,  to  search  out  the  total  number 
of  items  used  in  all  applications  for  a  particular  plane  or  set 


of  planes. 
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C •  Results  of  Statistical  Analysis 

First  the  pattern  of  obsolescence  was  studied.  By  pattern  of 
obsolescence,  we  mean  the  frequency  distributions  of  individual  items 
according  to  their  length  of  time-to-obsolescence  (or  life  span),  i.e., 
how  many  items  had  short,  medium,  or  long  life  span?  Thus,  we  are  in¬ 
terested  in  studying  the  pattern  of  obsolescence  that  characterizes 
groups  of  items.  Accordingly,  the  data  were  sorted  by  Federal  Stock 
Code  (FSC)  and  Technical  Supply  Maintenance  Code  (TSMC) • 

For  each  of  the  samples  thus  formed,  a  histogram  was  made  of  a 
distribution  of  life  spans.  Although  most  of  the  histograms  did  not 
exhibit  any  meaningful  pattern,  some  of  them  did  indicate  a  unimodal 
type  distribution.  Figures  1-4  show  a  series  of  histograms  which  are 
typical  of  this  latter  group. 

Those  probability  distributions  that  were  appropriate  to  this  type 

of  data  were  studied  in  order  to  find  if  any  could  be  fitted  to  the 

bt-afe*5t-U 

data.  It  was  found  that  the  failure  distribution  f(t)  ■  a  b  e  v  v 
(for  discussion  on  the  properties  of  this  distribution,  see  Chapter  V) 

closely  approximated  the  patterns  shown  by  the  histograms.  Trial  values 
of  the  parameters  a  and  b  were  tried  out  and  a  chi  square  test  perform¬ 
ed  to  indicate  whether  the  failure  distribution  fitted  the  distributions 
of  life  spans  for  particular  sets  of  parameters.  It  was  found  that 
the  distributions  of  life  spans  of  the  FSC-TSMC  combinations  1560-ADGA, 
1650-ADGA,  2810-PWAC  and  2840-PWAC  (Figures  1-4)  were  reasonably  de¬ 
scribed  by  failure  distributions.  The  parameters  a=l,  b«.0343  were  best 
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No.  of  Items  2  No.  of  Items 


Time-to-Obsolescence  (months) 

g.  1:  Frequency  Distribution  of  Time-to-Obsolescence  for  1560-ADGA  Sample* 
*FSN  1560:  Airframe  structural  components 
TSMC  ADGA:  Douglas  Aircraft  Company 


Time-to-Obsolescmce  (months) 

Fig.  2:  Frequency  Distribution  of  Time-to-Obsolescmce  for  1650-ADGA  Sample* 
*FSN  1650:  Aircraft  landing  gear  component 


No.  of  Items  ^  No.  of  Items 
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.  3:  Frequency  Distribution  of  Tlme-to-Obsolescence  for  2810-PWAC  Sample* 
*FSN  2810:  Gasoline  reciprocating  engines,  aircraft,  &  components 
TSMC  PWAL:  Wright  Aeronautical  Corporation 


Fig.  4:  Frequency  Distribution  of  Tlme-to-Obsolescence  for  2840-PWAC  Sample* 
*7SN  2840:  Gas  turbine  and  jet  engine  components 
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for  ADGA  distributions  and  a«l,  b-,1429  bast  for  PWAC  distributions.  Two 
particular  values  of  life  span,  116  and  118  from  ADGA  and  PWAC  distri¬ 
butions  respectively,  appeared  to  have  abnormally  high  frequencies,  so, 
were  ignored  when  fitting  distributions.  The  1650-ADGA  appeared  to  be 
short  of  its  actual  life  span,  but  it  resembled  a  truncated  1560-ADGA 
distribution  so  this  was  considered  when  fitting  distribution.  The 
2840- PWAC  distribution,  although  it  was  fitted,  had  frequencies  occur¬ 
ring  on  a  limited  number  of  values  of  life  span. 

Using  these  four  groups  of  data,  a  stepwise  regression  analysis  was 
carried  out  using  eight  quantitative  characteristics  of  each  obsolescent 
item  as  the  independent  variables  against  the  life  span  (y)  of  item  as 
dependent  variable.  The  eight  variables  are: 

*  unit  price 
Xj  “  lead  time 

*  total  Quantity  purchased 
X^  *  total  number  of  purchases 
X^  ■  number  of  applications 

Xg  «  number  of  units  per  application 
X7  *  maintenance  percent 
Xg  *  overhaul  percent 

Some  of  the  data  lacked  quantitative  information  for  lead  time, 
total  quantity  purchased  and  number  of  purchases.  These  data  were  not 
used  in  the  regression  analysis.  As  a  comparison,  a  regression  analysis 
was  also  done  on  ADGA  groups  where  frequencies  of  life  span  value  116  was 
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excluded.  Therefore,  a  total  of  six  stepwise  regression  analyses  were 
done.  Since  the  F  level  to  enter  a  variable  was  set  at  4.5,  no  in¬ 
dependent  variable  enters  the  step-wise  regression  analysis  in  two 
cases.  That  is  to  say,  none  of  the  eight  quantitative  variables  were 
found  to  be  statistically  related  to  items  life  span.  Results  of  the 
regression  analyses  are: 

1)  1560-ADGA  Group 

a.  Including  items  with  life  span  of  116  months. 

y  =  88.92314  +  0.00643X  -  3.66734X2  +  4.17063X3  +  2.95459X4 

(0.00272)  (1.14354)  (0.90705)  (1.29479) 

Standard  error  of  estimate  =  24.7458 
R  *  0.453 

b.  Excluding  items  with  life  span  of  116  months. 

y  *=  88.41669  +  0.00729X][  -  4.19640X,,  +  3.92748X3 
(0.00265)  (1.12390)  (0.92716) 

Standard  error  of  estimate  =  23.1396 
R  =  0,537 

2)  1650-ADGA  Group 

a.  Including  items  with  life  span  of  116  months, 
y  -  92.67528  +  5.65230X4  -  2.20044X5 

(1.26133)  (0.85662) 

Standard  error  of  estimate  »  15.4619 
R  -  0.687 
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b.  Excluding  items  with  life  span  of  116  months* 

No  independent  variable  entered  the  regression  analysis. 

3)  2810- PWAC  Group 

No  independent  variable  was  found  to  be  significant. 

4)  2840- PWAC  Group 

y  =  104.28982  -  1.27827X2  -  139.98844X3 
(0.47045)  (53.05405) 

Standard  error  of  estimate  =  14.2427 
R  -  0.466 

The  number  in  parenthesis  and  directly  below  each  regression  co¬ 
efficient  is  the  standard  deviation  of  the  estimate  of  coefficient. 

From  the  above  analysis,  the  following  tentative  conclusions  may 
be  drawn: 

1.  For  airframe  structural  components,  an  item  with  higher  unit 
value  tends  to  have  a  longer  life  span.  Perhaps  there  is  some  reluctance 
on  the  part  of  the  inventory  clerks  to  declare  these  high  valued  items 

as  obsolescent  items. 

2.  It  appears  that  the  longer  the  lead  time  the  more  likely  that 
obsolescence  will  set  in.  This  was  the  case  for  both  airframe  structural 
parts  and  engine  parts.  One  conjecture  is  that  the  long  lead  time  may 
mean  a  less  flexible  program.  Since  such  a  program  is  not  desirable, 
from  the  standpoint  of  an  efficient  supply  system,  there  is  a  tendency 

to  eliminate  the  item  concerned  from  the  system. 


3.  In  two  cases,  the  regression  analysis  indicates  that  the  total 
quantity  purchased  during  an  item's  life  span  is  positively  correlated 
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to  the  life  span  Itself,  However,  a  plausible  explanation  is  that  if 
an  item  has  a  long  life  span,  it  is  more  likely  that  more  purchases 
will  be  made. 

As  we  consider  how  declarations  of  obsolescence  have  been  made  here 
tofore  we  can  gain  some  appreciation  of  the  reasons  for  the  relatively 
weak  dependencies  of  time- to- obsolescence  upon  such  data  as  have  been 
available. 

There  are,  basically,  two  grounds  for  declaration  of  obsolescence 
of  line  items  in  Naval  inventory.  One  is  that  on  strategic  grounds  a 
weapons  system  is  declared  to  be  obsolete  after  some  point  in  time. 

That  is,  it  is  declared  that  from  some  date  on  this  system  will  be 
out  of  use.  The  other  is  that  a  stock  control  clerk  in  a  supply- 
demand  control  point,  following  some  rule  of  thumb  as  well  as  his  time 
permits,  and  his  assiduousness  compels,  notes  absence  of  demand  for 
items  over  a  period  of  time,  and  makes  a  declaration  of  obsolescence 
if  periods  of  sufficient  length  elapse  during  which  demand  is  sufficient 
ly  low  to  satisfy  the  rule  of  thumb. 

So  far  as  the  first  ground  for  declaration  of  obsolescence  is 
concerned,  this  is  less  a  basis  for  termination  of  the  system  as  a 
collection  of  line  items  in  Naval  inventory  than  a  basis  for  the  down¬ 
grading  of  the  system  in  importance  in  the  Navy's  arsenal  and  its  re¬ 
placement  by  some  other  system  or  systems.  This  decision  frequently 
commits  the  Navy  to  another  decision,  which  is  that  no  further  pro¬ 
curements  shall  be  made  in  support  of  the  older  system. 
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Concelvably,  the  application  of  models  of  the  sorts  developed  in  Chap 
ters  IV  &  VI  will  lead  to  a  reorganisation  of  the  structure  of  these 
decisions,  and  therefore,  to  the  entire  process  of  obsolescence 
generation. 

If  a  system  has  been  downgraded  it  will  continue  in  use  for  some 
period,  but  at  a  lower  use  rate.  This,  in  turn,  will  reduce  demand  for 
spares  items,  insofar  as  the  demand  is  due  to  wearout.  The  Bayesian 
feature  of  these  models  will  then  generate  higher  and  higher  probabili¬ 
ties  that  demand  is  running  at  lower  levels.  With  known  costs  of  disposal 
and  of  inventory  operation,  information  can  be  fed  to  high  level  echelons 
at  which  final  decisions  about  the  obsolescence  of  a  system  can  be  made 
definitively. 

When  the  matter  at  issue  is  the  simpler  case  of  a  decision  that  an 
item  can  be  declared  obsolete  by  a  stock  control  clerk  because  of  low 
demand,  this  entire  decision  can  be  handled  within  the  structure  of  the 
model,  on  a  more  thorough  and  reliable  basis.  Once  again,  the  costs  of 
carrying  the  item  on  inventory,  and  disposal  costs,  as  well  as  stockout 
costs,  would  be  taken  into  account  in  using  the  model. 

This  application  of  the  models  developed  in  this  study  to  generate 
obsolescence  decisions,  along  with  the  installation  of  appropriate  data 
processing  procedures  of  the  sort  sketched  out  in  an  accompanying  pro¬ 
posal,  would  lead  to  a  number  of  desirable  results:  sharper  definition 
of  obsolescence  determination,  less  confusion  in  record  keeping,  an  im¬ 


provement  in  the  quality  of  data  about  inventory  line  items.  As  a 
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consequence,  Chart  would  bo  developed  o  body  of  doe*  which  would  allow 
analysts  to  arrive  at  raliabla  and  neanintful  conclusions  as  to  tha 
dapandancias  of  tiae*to«obsoleaeeace  upon  such  data  as  characteristics 
of  the  itaa,  its  application,  and  tha  purchase  contracts  which  embody 
tha  orders. 


CHAPTER  VIII 
SUMMARY  AND  CONCLUSIONS 

This  chapter  consists  of  a  summary  of  the  work  which  has  been 
carried  out  under  this  contract,  and  some  recommendations  for  further 
work  which  are  spelled  out  in  detail  in  accompanying  proposals.  In  the 
course  of  this  work  an  attempt  was  made  to  unearth  meaningful  empirical 
relationships  between  obsolescence  and  engineering,  administrative  and 
program  data  relating  to  Naval  Aviation  parts  which  are  within  the  cog¬ 
nizance  of  the  U.S .  Naval  Aviation  Supply  Office.  In  addition,  a 
series  of  models  was  developed  and  explored  to  varying  degrees  which 
embody  more  and  more  sophisticated  notions  of  just  what  obsolescence  is 
operationally,  and  how  it  can  be  revealed  in  a  useful  fashion. 

It  has  been  found  that  by  virtue  of  the  nature  of  inventory  admin¬ 
istration  procedures  which  are  currently  in  use  at  ASO,  and  which  are 
probably  typical  of  procedures  carried  out  at  other  Supply-Demand-Con¬ 
trol  Points,  obsolescence  is  recorded  in  a  fashion  which  is  quite  un¬ 
predictable,  and  which  does  not  lend  itself  to  meaningful  statistical 
analysis.  Recording  of  obsolescence,  as  well  as  maintenance  of  records 
as  to  dates,  amounts  and  prices  of  purchases  of  items,  is  incomplete  and 
quite  dependent  on  such  variable  matters  as  work  loads,  assiduousness  of 
clerical  staff  and  their  supervisors,  availability  of  computing  and  com¬ 
puter  programming  facilities,  etc.  Moreover,  it  is  our  impression  that 
these  matters,  particularly  those  relating  to  computing  facilities,  are 
serious  bottlenecks. 

Finally,  even  if  these  unplanned  problems  were  cleared  away,  there 
remains  the  fundamental  difficulty  that  the  inventory  control  procedures 
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which  comprise  the  ideal,  the  plan  according  to  which  ASO  tries  to  oper¬ 
ate,  represent  a  pastiche  of  notions,  some  inconsistent  with  each  other, 
others  redundant,  with  no  real  attempt  to  face  a  set  of  system-wide  ob¬ 
jectives.  These  procedures  are  described  in  Appendix  C  to  this  report. 

A  comparison  of  this  appendix  with  Chapter  VIII  of  Whitin  [13]  ,  which 
was  published  in  1957,  is  most  illuminating.  The  chapter  referred  to 
contains  a  description  of  the  Navy  Inventory  control  system  which  was  in 
effect  prior  to  the  publication  of  that  book,  with  some  diagnosis  and 
recommendation  for  change. 

For  some  time  before  and  since  the  publication  of  Whitin's  book 
BuSandA  has  been  supporting  research  aimed  at  improvement  of  Navy  In¬ 
ventory  procedures.  However,  it  has  appeared  to  those  of  us  who  spent  a 
number  of  weeks  at  ASO  and  more  weeks  subsequently  digesting  what  we  had 
learned  there,  that  thus  far  the  effect  of  this  research  has  been  rela¬ 
tively  small  because  results  have  been  fed  into  existing  organizations 
in  small  amounts,  on  tentative  bases.  Recommended  procedures  have  been 
modified  so  substantially  in  order  to  accommodate  existing  organizational 
structures,  ways  of  doing  things  and  prejudices,  as  to  be  difficult  to 
recognize.  Moreover,  there  has  been  a  tendency  to  take  the  position  that 
what  is  good  in  one  environment  is  good  in  Isolation.  Hence,  a  collage 
of  measures,  techniques  and  procedures  has  been  assembled  from  widely 
scattered  sources  and  these  comprise  the  system. 

As  a  consequence'  ot  the  present  study  several  models  have  been  de¬ 
veloped  which  represent  successively  better  global  approaches  to  the 
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of  these  models  enables  the  skilled  Inventory  manager  to  take  ancillary 
Information,  or  what  the  Navy  refers  to  as  "program"  data,  into  account 
in  a  systematic  fashion.  Specifically,  if  on  budgetary  or  strategic 
grounds  a  particular  weapons  system  is  expected  to  become  less  signifi¬ 
cant  in  the  activities  of  the  Atlantic  Fleet,  and  an  aircraft  type  will 
be  flying  fewer  hours,  then  a  priori  distributions  of  the  probability  of 
obsolescence,  or  a  priori  probabilities  of  being  in  one  or  another  state 
of  demand  can  be  altered  to  reflect  that  expectation.  Moreover,  these 
models  are  flexible  as  to  lead  time,  and  other  inventory  parameters. 

It  is  true  that  these  models  make  certain  rather  stringent  demands 
of  the  inventory  manager  which  cannot  be  solved  simply  by  improving  the 
quality  of  data  handling.  All  of  these  demands  are  in  the  area  of  ac¬ 
quiring  better  understanding  of  the  cost  nature  of  various  conditions  of 
the  inventory.  How  much  does  it  really  cost  the  Navy  to  be  out  of  stock 
of  a  particular  part?  How  much  does  it  really  cost  to  carry  such  and 
such  an  item  in  inventory,  to  order  it,  etc.?  These  questions  should  be 
answered,  somehow  or  other.  However,  even  if  the  first-cut  answers  are 
semi-educated  guesses,  they  will  enable  the  models  to  be  operated,  and 
in  fact  these  answers  could  be  checked  in  terms  of  the  resultant  in¬ 
ventory  policies  which  would  be  generated. 

What  is  suggested,  then,  is  that  these  models  can  eventually  serve 
as  inventory  control  procedures.  Indeed,  proposals  for  further  study 
and  for  installation  accompany  this  report.  However,  what  is  perhaps 
more  significant  immediately  is  that  these  models  could  serve  as 

_  _ 
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measuring  instruments  for  the  development  of  measures  of  some  of  the 
more  fugitive  inventory  parameters  of  which  theoreticians  speak  so 
casually,  such  as  stockout  cost,  setup  cost,  etc.  Thus,  a  sensible 
procedure  which  might  be  adopted  could  involve  four  phases  (discussed 
in  detail  in  the  accompanying  proposal  for  further  work): 

1.  Set  up  and  begin  operation  of  a  computer  system  designed 

to  look  up  appropriate  control  policies  for  separate  items. 

These  policies  would  be  generated  elsewhere  according  to 
some  model. 

2.  At  the  same  time,  continue  exploration  of  models.  Eventually, 
tables  should  b-;-  generated  from  the  appropriate  model  or 
models  which  would  be  available  for  lookup  in  inventory  control. 

3.  Load  into  the  machine  system,  as  they  are  developed,  tables 
based  on  successively  more  satisfactory  models.  Running 
parallel  to  the  existing  inventory  control  system,  run  this 
new  system  and  compare  results.  As  confidence  in  these  tables 
developes,  studies  are  begun,  using  these  models,  which  would 
lead  to  measures  of  inventory  system  parameters  consistent  with 
current  policies,  (i.e.  policies  in  use  before  any  of  these 
changes  have  begun).  This  procedure,  iterated  and  played 
against  inventory  managers  and  those  responsible  for  high 
level  logistical  policy,  would  lead,  eventually,  to  accept¬ 
able  measures  of  such  parameters  as  setup  cost,  lead  time, 
stockout  cost,  etc.  In  addition,  during  this  phase,  studies 
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could  be  carried  out  which  would  reveal,  in  detail,  the 
amount  of  savings  in  physical  inventory  costs  which  would 
accrue  from  the  installation  of  the  system,  as  well  as  the 
improvement  in  quality  of  inventory  data.  Note  that  up  to 
this  time  the  normal  operations  of  the  SDCP  have  not  been 
interfered  with.  All  the  work  described  in  these  three 
phases  has  been  carried  on  outside  the  current  operation. 

4.  Once  satisfaction  has  been  achieved  with  respect  to  the 

utility  of  the  new  system,  it  could  be  phased  in  gradually, 
while  the  old  one  is  phased  out.  Conceivably,  with  low 
demand,  high  value  items,  or  with  items  with  very  high 
stockout  costs  or  other  special  problems,  phase-in  might 
be  a  much  slower  process.  Indeed,  it  might  be  desirable 
that  some  items  never  be  completely  controlled  by 
mechanized  procedures,  but  be  constantly  subject  to  review 
by  highly  skilled  inventory  specialists.  However,  for  the 
bulk  of  the  400,000  items  within  the  cognizance  of  ASO,  as 
well,  probably,  as  the  remainder  of  the  1.2  million  items 
handled  by  other  SDCP's  for  the  Navy,  this  sort  of  proce¬ 
dure  would  probably  turn  out  to  be  a  most  economical  and 
satisfactory  one. 

Although  it  will  not  be  gone  into  detail  here  conservative  prelimi¬ 
nary  estimates  indicate  that  inventory  administration  costs  (i.e.  costs 


of  maintaining  files,  ordering,  disposing,  etc.)  for  the  400,000  items 
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handled  by  ASO,  could  probably  be  reduced  by  a  substantial  amount  each 
year.  This  does  not  take  into  account  the  actual  physical  inventory 
savings  due  to  changes  in  amounts  of  inventory  held,  changes  in  numbers 
of  stockouts,  disposals,  orders,  etc.  What  is  being  asserted  here  is 
that  the  costs  of  placing  exactly  the  same  orders,  carrying  exactly  the 
same  amounts  of  all  items,  disposing  of  the  same  amounts,  etc.  at  the 
same  times,  would  be  much  less  than  they  now  are.  Savings  resulting 
from  changing  these  policies  would,  in  all  likelihood,  be  much  more  sub¬ 
stantial,  but  their  magnitude  would  begin  to  be  revealed  only  during  the 
third  phase  of  the  four  phase  procedure  laid  out  above. 


APPENDIX  A 


POLICIES  RESULTING  FROM  DYNAMIC  PROCUREMENT  - 
DISPOSAL  MODEL  RUNS 

Before  presenting  the  results  of  the  computer  runs,  we  shall  first 
describe  how  the  values  of  parameters,  which  characterize  the  demand  and 
obsolescence  distributions  and  cost  functions,  are  varied  in  each  of  the 
runs . 

Each  run  is  identified  by  a  three-digit  number.  The  first  digit 
identifies  one  of  the  five  different  sets  of  assumptions  regarding  the 
demand  and  obsolescence  distributions.  The  last  two  digits  identify  a 
different  cost  combination.  These  identification  codes  are  indicated  in 
Tables  1-A  and  2-A. 

Policies  resulting  from  each  computer  run  are  tabulated  in  the 
remainder  of  this  appendix.  Note  that  for  each  run,  we  assumed  an  in¬ 
ventory  process  with  a  finite  horizon  consisting  of  10  periods. 


TABLE  1-A 


Variations  in  Demand  and 
Obsolescence  Distributions 


Run  # _ Demand  Distribution _ Obsolescence  Distribution 


1XX 

Poisson  demand  with 
mean  of  30;  demand 
rate  not  changing 
through  time 

Exponential  distribution 
with  obsolescence  rate  of 
10%. 

2XX 

Poisson  demand  with 
decreasing  mean  de¬ 
mand  rate  through 
time  as  follows: 

Same  as  in  1XX  runs . 

Time  Mean  demand 

1,2  40 

3,4  35 

5,6  30 

7,8  25 

9,10  20 

3XX 

Uniform  demand  with 
mean  of  30;  demand 
rate  not  changing 
through  time 

Same  as  in  1XX  runs . 

4XX 

Same  as  in  3XX  runs 

Uniform  distributions: 

10%  obsolescence  rate. 

5XX 

Uniform  distribution 
decreasing  mean  de¬ 
mand  rate  as  in  2XX 

Same  as  in  4XX  runs . 

runs 
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TABLE  2 -A 

Variations  in  Cost  Parameters 


Run  # 

*  ** 

Disposal 

Cost 

•k 

Fixed  Ordering 
Cost 

*  ** 

Storage 

Cost 

* 

Out-of-stock 

Cost 

X  01 

0-2 

Z 

10 

0.1  Z 

5  Z 

X  02 

0-2 

Z 

10 

0.1  Z 

30  Z 

X  03 

0-2 

Z 

40 

0.05  Z 

10  Z 

X  04 

0-2 

Z 

40 

0.05  Z 

30  Z 

X  05 

0-2 

Z 

40 

0.1  Z 

30  Z 

X  06 

15-2 

Z 

10 

0.1  Z 

5  Z 

X  07 

15-2 

Z 

10 

0.2  Z 

5  Z 

X  08 

15-2 

Z 

10 

0.1  Z 

30  Z 

X  09 

15-2 

z 

40 

0.1  Z 

5  Z 

X  10 

15-2 

z 

40 

0.1  Z 

30  Z 

X  11 

0-2 

z 

40 

0.02  Z 

10  Z 

X  12 

0-2 

z 

40 

0.02  Z 

30  Z 

* 

All 

costs 

are  expressed  as  multiples 

of  unit  price 

of  the  stock. 

(i.e.,  unit  price  ■  1) 

Z  =  amount  to  be  disposed,  stored  or  to  which  a  penalty  cost  is 
charged . 
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TABLE  3 -A 

Policies  Resulting  from  Model  Runs 


Run  101 


** 

Period 

h 

f2 

Ik 

10 

27 

34 

43 

44 

9 

29 

3L 

77 

78 

8 

29 

38 

108 

109 

7 

29 

38 

137 

138 

6 

29 

38 

165 

166 

5 

29 

38 

192 

193 

4 

29 

38 

216 

217 

3 

29 

38 

238 

239 

2 

29 

38 

256 

257 

1 

29 

38 

261 

262 

Run  103 

•ifir 

Period 

_h 

fl 

h 

\ 

10 

27 

37 

47 

48 

9 

29 

68 

82 

83 

8 

27 

70 

115 

116 

7 

28 

70 

146 

147 

6 

27 

70 

176 

177 

5 

28 

70 

206 

207 

4 

28 

70 

235 

236 

3 

28 

70 

263 

264 

2 

28 

70 

290 

291 

1 

28 

70 

. 

. 

Run  102 


irk 

Period 

Jk 

h 

h 

h 

10 

34 

39 

46 

47 

9 

35 

42 

81 

82 

8 

35 

42 

112 

113 

7 

35 

42 

142 

143 

6 

35 

42 

170 

171 

5 

35 

42 

196 

197 

4 

35 

42 

221 

222 

3 

35 

42 

242 

243 

2 

35 

42 

260 

261 

1 

35 

42 

265 

266 

Run  104 

Period** 

h 

h 

h 

h 

10 

31 

39 

47 

48 

9 

33 

70 

84 

85 

8 

32 

73 

117 

118 

7 

32 

74 

148 

149 

6 

32 

73 

179 

180 

5 

32 

74 

208 

209 

4 

32 

73 

235 

236 

3 

32 

74 

265 

266 

2 

32 

73 

292 

293 

1 

32 

74 

- 

- 

* 

Inventory  Level.  H _ Policy 

H  <  Hx  Buy  Kx  -  H 

H  >  Dispose  H  - 

** 


Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n. 
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TABLE  3 -A  Cont'd 

Policies  Resulting  from  Model  Runs 
Run  105  Run  106 


** 

Period 

111 

i 

h 

h 

ikic 

Period 

i 

h 

10 

31 

39 

46 

47 

10 

25 

32 

41 

199 

9 

32 

70 

82 

83 

9 

28 

37 

73 

164 

8 

32 

73 

114 

115 

8 

28 

37 

104 

173 

7 

32 

73 

145 

146 

7 

28 

37 

133 

192 

6 

32 

73 

174 

175 

6 

28 

37 

161 

215 

5 

32 

73 

202 

203 

5 

28 

37 

188 

240 

4 

32 

73 

229 

230 

4 

28 

37 

213 

266 

3 

32 

73 

254 

255 

3 

28 

37 

236 

290 

2 

32 

73 

276 

277 

2 

28 

37 

- 

- 

* 

1 

32  (?) 

73(?)  *** 

. »  . »  . » . 

Tnf  x 

1 

28 

37 

- 

- 

Run  107 
** 

Period 

i 

i 

h 

i 

Run  108 

Period 

!i 

fl 

h 

10 

25 

32 

39 

118 

10 

33 

38 

45 

203 

9 

28 

37 

71 

117 

9 

35 

42 

79 

169 

8 

28 

37 

100 

137 

8 

35 

42 

110 

178 

7 

28 

37 

127 

161 

7 

35 

42 

139 

198 

6 

28 

37 

153 

186 

6 

35 

42 

167 

221 

5 

28 

37 

174 

211 

5 

35 

42 

194 

246 

4 

28 

37 

189 

233 

4 

35 

42 

219 

271 

3 

28 

37 

189 

250 

3 

35 

42 

242 

296 

2 

28 

37 

189 

258 

2 

35 

42 

- 

- 

1 

28 

37 

189 

259 

1 

35 

42 

- 

- 

* 

Inventory  Level. 

H 

Policy 

H<  Hx 

Buy 

-  H 

H>  H 

Dispose  H  - 

** 

Real  time  periods; 

i.e.  period  n  + 

1  is  later 

than  period 

n. 

*** 

Through 

an  operator  error 

,  the  full  output  for  this  period  was 

not 

printed.  and  values  were  inferred  from  the  partial  output,  but 
the  and  K„  values  could  not  be  so  inferred. 


TABLE  3 -A  Cont'd 

Policies  Resulting  from  Model  Runs 


Run  109 


** 

Period 

h 

h 

f2 

^2 

10 

23 

35 

44 

202 

9 

25 

65 

77 

167 

8 

22 

67 

108 

176 

7 

24 

67 

138 

196 

6 

23 

67 

167 

220 

5 

24 

67 

195 

246 

4 

23 

67 

222 

272 

3 

24 

67 

248 

298 

2 

23 

67 

- 

- 

1 

24 

67 

- 

- 

Run  111 
** 

Period 

i 

h 

*2 

^2 

10 

27 

38 

48 

49 

9 

29 

68 

85 

86 

8 

27 

95 

119 

120 

7 

28 

71 

151 

152 

6 

28 

71 

183 

184 

5 

28 

71 

214 

215 

4 

28 

71 

244 

245 

3 

28 

71 

273 

274 

2 

28 

71 

- 

- 

1 

28 

71 

- 

- 

Inventory  Level, H  Policy 


H  <  Hx 

Buy  Kx  -  H 

H  >  H2 

Dispose  H  - 

Real  time  periods;  i.e.  period  n 


Run  110 


irk 

Period 

i 

h 

i 

i 

10 

31 

39 

46 

203 

9 

32 

70 

80 

170 

8 

32 

73 

112 

180 

7 

32 

73 

142 

200 

6 

32 

73 

171 

224 

5 

32 

73 

199 

250 

3 

32 

73 

227 

277 

3 

32 

73 

- 

- 

2 

32 

73 

- 

- 

1 

32 

73 

- 

- 

Run  112 

irk 

Period 

h 

f2 

10 

31 

39 

48 

49 

9 

33 

71 

86 

87 

8 

32 

97 

121 

122 

7 

33 

74 

153 

154 

6 

32 

74 

185 

186 

5 

32 

74 

216 

217 

4 

32 

74 

246 

247 

3 

32 

74 

275 

276 

2 

32 

74 

- 

- 

1 

32 

74 

1  is  later 

than 

period 

n. 
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TABLE  3 -A  Cont'd 

Policies  Resulting  from  Model  Runs 
Run  201  Run  202 


•Me 

Period 

H1 

K1 

K2 

H2 

** 

Period 

H1 

K1 

h 

h 

10 

17 

23 

31 

32 

10 

23 

27 

33 

34 

9 

18 

26 

54 

55 

9 

24 

30 

57 

58 

8 

23 

32 

81 

82 

8 

29 

36 

84 

85 

7 

23 

32 

105 

106 

7 

29 

36 

109 

110 

6 

29 

38 

134 

135 

6 

35 

42 

138 

139 

5 

29 

38 

162 

163 

5 

35 

42 

165 

166 

4 

34 

43 

193 

194 

4 

41 

48 

196 

197 

3 

34 

43 

222 

223 

3 

41 

48 

225 

226 

2 

39 

49 

253 

254 

2 

46 

54 

257 

258 

1 

39 

49 

281 

282 

1 

46 

54 

284 

285 

Run  203 
** 

Period 

H1 

K1 

CM 

5* 

«2 

Run  204 
Period** 

H1 

i 

h 

i 

10 

17 

26 

34 

35 

10 

20 

28 

35 

36 

9 

19 

47 

58 

59 

9 

22 

49 

60 

61 

8 

23 

70 

87 

88 

8 

27 

72 

88 

89 

7 

23 

61 

113 

114 

7 

27 

64 

115 

116 

6 

28 

65 

145 

146 

6 

32 

68 

146 

147 

5 

28 

90 

175 

176 

5 

32 

74 

177 

178 

4 

33 

76 

209 

210 

4 

38 

79 

211 

212 

3 

33 

81 

243 

244 

3 

37 

85 

245 

246 

2 

38 

86 

281 

282 

2 

43 

90 

282 

283 

1 

37 

92 

- 

1 

43 

95 

“ 

- 

* 

Inventory  Level, 

H 

Policy 

H  <  H  Buy  Kx  -  H 

H>  H2  Dispose  H  -  K2 

Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n. 


I 

I 
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TABLE  3 -A  Cont'd 

Policies  Resulting  from  Model  Runs 


Run  205 

irk 

Period 

i 

h 

h 

i 

» - - 

Run  206 

•kie 

Period 

h 

h 

S 

10 

20 

28 

34 

35 

10 

15 

21 

28 

185 

9 

22 

49 

58 

59 

9 

18 

26 

51 

140 

8 

27 

71 

86 

87 

8 

23 

32 

77 

145 

7 

27 

63 

112 

113 

7 

23 

31 

101 

160 

6 

32 

68 

142 

143 

6 

28 

37 

130 

183 

5 

32 

74 

171 

172 

5 

28 

37 

157 

209 

4 

38 

78 

204 

205 

4 

33 

43 

189 

240 

3 

37 

84 

235 

236 

3 

33 

43 

219 

271 

2 

43 

89 

269 

270 

2 

38 

48 

- 

- 

1 

42 

95 

- 

- 

1 

38 

48 

- 

- 

Run  207 
** 

Period 

H1 

K1 

K2 

h 

105 

Run  208 
Period** 

Hi 

22 

fl 

27 

h 

^2 

10 

15 

21 

27 

10 

32 

189 

9 

17 

25 

49 

94 

9 

24 

30 

56 

145 

8 

23 

31 

73 

109 

8 

29 

36 

82 

150 

7 

23 

31 

96 

127 

7 

29 

36 

106 

165 

6 

28 

37 

123 

155 

6 

35 

42 

135 

188 

5 

28 

37 

147 

181 

5 

35 

42 

162 

214 

4 

33 

42 

173 

210 

4 

41 

48 

194 

245 

3 

33 

42 

192 

238 

3 

41 

48 

223 

276 

2 

38 

48 

206 

265 

2 

46 

54 

- 

- 

1 

38 

48 

221 

285 

1 

46 

54 

- 

- 

Inventory  Level.  H  Policy _ 

H  <  Hx  Buy  Kx  -  H 

H  >  Dispose  H  -  Kj 

** 


Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n 
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TABLE  3 -A  Cont'd 

•jf 

Policies  Resulting  from  Model  Runs 


Run  209 


** 


Period 

Hi 

K1 

K2 

H2 

10 

13 

24 

32 

189 

9 

16 

44 

54 

143 

8 

19 

66 

80 

148 

7 

18 

58 

106 

164 

6 

24 

62 

135 

188 

5 

23 

68 

164 

214 

4 

29 

72 

197 

246 

3 

28 

78 

228 

278 

2 

33 

83 

- 

- 

1 

32 

88 

- 

- 

Run  211 

Period 

H1 

K1 

K2 

»2 

10 

17 

26 

35 

36 

9 

19 

47 

61 

62 

8 

23 

70 

90 

91 

7 

23 

61 

118 

119 

6 

28 

66 

150 

151 

5 

28 

91 

181 

182 

4 

33 

76 

217 

218 

3 

33 

82 

252 

253 

2 

38 

87 

292 

293 

1 

38 

92 

* 

• 

* 

Inventory  Level, 

H 

Policy 

Run  210 


ific 

Period 

i 

h 

h 

H2 

10 

20 

27 

33 

190 

9 

22 

48 

57 

145 

8 

27 

71 

84 

151 

7 

26 

63 

109 

167 

6 

32 

68 

139 

192 

5 

32 

74 

168 

218 

4 

38 

78 

200 

250 

3 

37 

84 

232 

282 

2 

43 

89 

- 

- 

1 

42 

95 

- 

- 

Run  212 
** 

Period 

Hi 

h 

i 

h. 

10 

20 

28 

35 

36 

9 

22 

49 

62 

63 

8 

27 

72 

91 

92 

7 

27 

64 

119 

120 

6 

33 

68 

152 

153 

5 

32 

95 

183 

184 

4 

38 

79 

219 

220 

3 

37 

85 

254 

255 

2 

43 

90 

294 

295 

1 

43 

96 

- 

. 

period  n. 


H  <  Hx 

H  >  h2 


Buy  Kx  -  H 
Dispose  H  - 


** 


Real  time  periods;  i.e.  period  n  +  1  is  later  than 
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TABLE  3- A  Cont'd 


Policies 

Resulting  from  Model  Runs* 

Run  301 

Run  302 

Period** 

h 

i 

Period** 

h 

h 

h 

10 

30 

44 

60 

61 

10 

50 

57 

60 

61 

9 

41 

56 

110 

111 

9 

53 

60 

116 

117 

8 

41 

57 

147 

148 

8 

53. 

60 

156 

157 

7 

41 

57 

174 

175 

7 

53 

60 

182 

183 

6 

41 

57 

199 

200 

6 

53 

60 

206 

207 

5  ^ 

41 

57 

220 

221 

5 

53 

60 

226 

227 

"4 

41 

57 

237 

238 

4 

53 

60 

243 

244 

3 

41 

57 

251 

252 

3 

53 

60 

257 

258 

2 

41 

57 

262 

263 

2 

53 

60 

267 

268 

1 

41 

57 

270 

271 

1 

53 

60 

274 

275 

Run  303 

Run  304 

Period** 

i 

^2 

Period** 

i 

fl 

h 

i 

10 

33 

53 

60 

61 

10 

45 

58 

60 

61 

9 

39 

60 

117 

118 

9 

48 

62 

118 

119 

8 

39 

81 

162 

163 

8 

48 

82 

167 

168 

7 

39 

101 

199 

200 

7 

48 

109 

204 

205 

6 

39 

102 

228 

229 

6 

48 

110 

233 

234 

5 

39 

102 

256 

257 

5 

48 

110 

260 

261 

4 

39 

102 

281 

282 

4 

48 

110 

285 

286 

3 

39 

102 

- 

- 

3 

48 

110 

- 

- 

2 

39 

102 

- 

- 

2 

48 

110 

- 

- 

1 

39 

102 

- 

- 

1 

48 

110 

- 

- 

•ff 

Inventory  Level.  H _ Policy 

H<  Ht  Buy  -  H 

H>  H2  Dispose  H  -  K2 

**Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n. 
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TABLE  3- A  Cont'd 

Policies  Resulting  from  Model  Runs* 


Run  305 


Period** 

i 

h 

s 

10 

45 

58 

60 

61 

9 

47 

60 

116 

117 

8 

47 

71 

162 

163 

7 

47 

94 

191 

192 

6 

47 

109 

216 

217 

5 

47 

109 

239 

240 

4 

47 

109 

260 

261 

3 

47 

109 

277 

278 

2 

47 

109 

291 

292 

1 

47 

109 

- 

- 

Run  307 

Period** 

h 

h 

^2 

10 

28 

41 

58 

133 

9 

39 

54 

99 

149 

8 

40 

55 

124 

175 

7 

40 

55 

147 

198 

6 

40 

55 

165 

219 

5 

40 

55 

181 

236 

4 

40 

55 

192 

250 

3 

40 

55 

198 

261 

2 

40 

55 

201 

267 

1 

40 

55 

202 

271 

* k 

Inventory  Level 

i  H 

Policy 

H  <  Hx 

Buy  K1 

-  H 

H  >  Hj  Dispose  H 


Run  306 


Period** 

h 

h 

h 

10 

28 

42 

59 

207 

9 

40 

55 

105 

194 

8 

40 

56 

137 

213 

7 

40 

56 

164 

236 

6 

40 

56 

189 

260 

5 

40 

56 

210 

283 

4 

40 

56 

- 

- 

3 

40 

56 

- 

- 

2 

40 

56 

- 

- 

1 

40 

56 

- 

- 

Run  308 

Period** 

h 

h 

10 

50 

56 

60 

208 

9 

53 

60 

114 

200 

8 

53 

60 

146 

222 

7 

53 

60 

172 

244 

6 

53 

60 

196 

268 

5 

53 

60 

218 

291 

4 

53 

60 

- 

- 

3 

53 

60 

- 

- 

2 

53 

60 

- 

- 

1 

53 

60 

_ 

**  Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n 
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TABLE  3- A  Cont’d 

Policies  Resulting  from  Model  Runs* 


Run  309  Run  310 


Period** 

h 

h 

i 

i 

Period** 

1 

h 

i 

10 

19 

46 

60 

208 

10 

44 

57 

60 

208 

9 

31 

60 

109 

197 

9 

47 

60 

114 

200 

8 

29 

67 

143 

217 

8 

47 

71 

153 

228 

7 

30 

89 

171 

242 

7 

47 

93 

180 

251 

6 

30 

91 

197 

268 

6 

47 

109 

207 

277 

5 

30 

91 

221 

293 

5 

47 

109 

- 

- 

4 

30 

91 

- 

- 

4 

47 

109 

- 

- 

3 

30 

91 

- 

- 

3 

47 

109 

- 

- 

2 

30 

91 

- 

- 

2 

47 

109 

- 

- 

1 

30 

91 

- 

- 

1 

47 

109 

- 

- 

Run  311 

Run  312 

Period** 

i 

h 

Period** 

i 

fl 

^2 

10 

32 

52 

- 

- 

10 

45 

57 

- 

- 

9 

39 

60 

- 

- 

9 

48 

62 

- 

- 

8 

39 

85 

- 

- 

8 

48 

87 

- 

- 

7 

39 

102 

- 

- 

7 

48 

110 

- 

- 

6 

39 

103 

- 

- 

6 

48 

111 

- 

- 

5 

39 

103 

- 

- 

5 

48 

110 

- 

- 

4 

39 

103 

- 

- 

4 

48 

110 

- 

- 

3 

39 

103 

- 

- 

3 

48 

110 

- 

- 

2 

39 

103 

- 

- 

2 

48 

110 

- 

- 

1 

39 

103 

- 

- 

1 

48 

110 

- 

- 

* 

Inventory  Level, 

H 

Policy 

H  <  Hx 

H  >  H2 

Buy  K  -  H 

Dispose  H  - 

K2 

**  Real 

time  periods; 

i.e. 

period 

n  +  1  is  later 

than  period 

n. 
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TABLE  3- A  Cont'd 


Policies  Resulting  from  Model  Runs* 
Run  401  Run  402 


Period** 

h 

h 

h 

Period** 

h 

h 

10 

30 

44 

60 

61 

10 

50 

57 

60 

61 

9 

36 

51 

107 

108 

9 

52 

59 

114 

115 

8 

39 

54 

137 

138 

8 

52 

59 

145 

146 

7 

40 

56 

160 

161 

7 

53 

60 

167 

168 

6 

41 

57 

182 

183 

6 

53 

60 

189 

190 

5 

42 

57 

204 

205 

5 

53 

60 

209 

210 

4 

42 

58 

223 

224 

4 

53 

60 

229 

230 

3 

42 

58 

243 

244 

3 

53 

60 

248 

249 

2 

43 

59 

261 

262 

2 

53 

60 

265 

266 

1 

43 

59 

278 

279 

1 

53 

60 

282 

283 

Run  403 

Run  404 

Period** 

\ 

h 

^2 

h 

Period** 

\ 

h 

K, 

10 

33 

53 

60 

61 

10 

45 

58 

60 

61 

9 

37 

58 

115 

116 

9 

47 

60 

117 

118 

8 

38 

60 

157 

158 

8 

47 

60 

162 

163 

7 

38 

60 

187 

188 

7 

47 

60 

191 

192 

6 

39 

96 

214 

215 

6 

48 

97 

218 

219 

5 

41 

104 

241 

242 

5 

49 

111 

245 

246 

4 

41 

107 

266 

267 

4 

49 

113 

270 

271 

3 

42 

110 

291 

292 

3 

49 

114 

295 

296 

2 

42 

112 

- 

- 

2 

50 

116 

- 

- 

1 

42 

114 

- 

- 

1 

50 

117 

- 

- 

•fa 

Inventory  Level.  H _ Policy 

H  <  Hx  Buy  Kx  -  H 

H  >  H2  Dispose  H  -  K2 

**  Real  time  periods;  i.e,  period  n  +  1  is  later  than  period  n. 
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I 

I 

I 

I 


TABLE  3-A  Cont'd 


Policies  Resulting  from  Model  Runs* 
Run  405  Run  406 


Period** 

1 

h 

h 

H2 

Period** 

i 

fl 

h 

10 

45 

58 

60 

61 

10 

28 

42 

59 

207 

9 

46 

60 

115 

116 

9 

35 

50 

102 

209 

8 

47 

60 

153 

154 

8 

38 

53 

129 

219 

7 

47 

60 

174 

175 

7 

40 

55 

153 

233 

6 

47 

73 

200 

201 

6 

40 

56 

175 

251 

5 

48 

109 

222 

223 

5 

41 

57 

197 

270 

4 

49 

112 

244 

245 

4 

42 

58 

217 

289 

3 

49 

113 

265 

266 

3 

42 

58 

- 

- 

2 

49 

114 

285 

286 

2 

42 

58 

- 

_ 

1 

49 

115 

- 

- 

1 

42 

59 

- 

- 

Run  407 

Run  408 

Period** 

i 

h 

h 

i 

Period** 

i 

h 

h 

10 

28 

41 

58 

133 

10 

50 

56 

60 

208 

9 

34 

49 

95 

155 

9 

52 

58 

112 

216 

8 

37 

52 

113 

171 

8 

52 

59 

137 

228 

7 

39 

54 

135 

190 

7 

53 

60 

161 

241 

6 

40 

55 

153 

207 

6 

53 

60 

182 

259 

5 

40 

56 

168 

223 

5 

53 

60 

204 

278 

4 

41 

57 

183 

238 

4 

53 

60 

223 

296 

3 

41 

57 

197 

252 

3 

53 

60 

- 

- 

2 

41 

57 

209 

265 

2 

53 

60 

- 

- 

1 

42 

58 

220 

277 

1 

53 

60 

Inventory  Level.  H 


Policy 


H  <  Hx  Buy  Kx  -  H 

H  >  H2  Dispose  H  -  ^ 


**  Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n 
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TABLE  3-A  Cont'd 

Policies  Resulting  from  Model  Runs* 
Run  409  Run  410 


Period** 

1 

h 

h 

Period** 

1 

h 

h 

h 

10 

19 

46 

60 

208 

10 

44 

57 

60 

208 

9 

26 

54 

106 

212 

9 

46 

59 

113 

217 

8 

28 

58 

136 

224 

8 

47 

60 

145 

234 

7 

29 

60 

160 

239 

7 

47 

60 

169 

248 

6 

30 

72 

184 

259 

6 

47 

72 

193 

268 

5 

32 

94 

207 

279 

5 

48 

109 

216 

288 

4 

32 

98 

229 

300 

4 

49 

112 

«• 

- 

3 

33 

101 

- 

- 

3 

49 

113 

m 

- 

2 

33 

104 

- 

- 

2 

49 

114 

m 

- 

1 

34 

106 

- 

- 

1 

49 

115 

- 

- 

Run  411 

Run  412 

Period** 

1 

h 

^2 

Period** 

h. 

h 

CM 

10 

32 

52 

10 

45 

57 

«• 

9 

36 

58 

9 

46 

60 

m 

8 

37 

60 

8 

47 

60 

m 

7 

38 

60 

7 

47 

60 

m 

6 

39 

100 

6 

48 

108 

«■ 

5 

41 

105 

5 

49 

112 

4 

41 

109 

4 

49 

114 

3 

42 

111 

- 

3 

49 

115 

2 

42 

114 

- 

2 

50 

117 

• 

- 

1 

43 

116 

• 

- 

1 

50 

119 

. 

Inventory  Level.  H  Policy 
H  <  Hx  Buy  IL  -  H 

H  >  H2  Dispose  H  -  K2 

**  Reel  tlmti  periods;  l.e.  period  n  +  1  Is  later  then  period  n. 
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TABLE  3-A  Cont'd 

Policies  Resulting  from  Model  Runs* 
Run  501  Run  502 


Period** 

i 

h 

h 

h 

Period** 

fl 

h 

h 

h 

10 

19 

30 

40 

41 

10 

33 

38 

40 

41 

9 

23 

35 

72 

73 

9 

34 

39 

77 

78 

8 

31 

45 

100 

101 

8 

43 

49 

106 

107 

7 

33 

47 

123 

124 

7 

44 

50 

128 

129 

6 

41 

56 

151 

152 

6 

53 

60 

155 

156 

5 

42 

57 

177 

178 

5 

53 

60 

181 

182 

4 

50 

67 

208 

209 

4 

63 

70 

213 

214 

3 

50 

68 

237 

238 

3 

63 

70 

241 

242 

2 

59 

77 

271 

272 

2 

72 

80 

275 

276 

1 

59 

78 

- 

- 

1 

72 

80 

- 

- 

Run  503 

Run  504 

Period** 

1 

h 

s 

Period** 

i 

h 

f2 

^2 

10 

20 

37 

40 

41 

10 

29 

39 

40 

41 

9 

22 

40 

78 

79 

9 

29 

40 

79 

80 

8 

31 

50 

114 

115 

8 

39 

50 

118 

119 

7 

31 

57 

142 

143 

7 

39 

57 

145 

146 

6 

40 

88 

176 

177 

6 

48 

91 

179 

180 

5 

41 

103 

206 

207 

5 

49 

110 

209 

210 

4 

49 

115 

244 

245 

4 

58 

121 

247 

248 

3 

50 

125 

278 

279 

3 

58 

131 

282 

283 

2 

58 

136 

- 

- 

2 

67 

142 

- 

- 

1 

59 

146 

- 

- 

1 

68 

152 

- 

- 

Inventory  Level.  H 


Police 


H  <  Hx  Buy  -  H 

H  >  H2  Dispose  H  -  K2 


**  Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n. 


A-17 


TABLE  3-A  Cont'd 

Policies  Resulting  from  Model  Runs* 


Run  505 

Period** 

h 

10 

28 

39 

9 

29 

40 

8 

38 

50 

7 

38 

50 

6 

48 

71 

5 

48 

109 

4 

57 

120 

3 

58 

130 

2 

67 

140 

1 

67 

150 

Run  507 

Period** 

h 

fl 

10 

16 

27 

9 

21 

33 

8 

30 

43 

7 

31 

45 

6 

40 

55 

5 

40 

56 

4 

49 

65 

3 

49 

66 

2 

58 

75 

1 

58 

76 

Run  506 

h 

i 

Period** 

40 

41 

10 

77 

78 

9 

113 

114 

8 

134 

135 

7 

165 

166 

6 

194 

195 

5 

227 

228 

4 

257 

258 

3 

294 

295 

2 

- 

- 

1 

Run  508 

h 

h 

Period** 

39 

113 

10 

63 

120 

9 

83 

135 

8 

102 

152 

7 

127 

177 

6 

147 

198 

5 

172 

226 

4 

194 

250 

3 

220 

279 

2 

. 

_ 

1 

i 

h 

16 

28 

39 

186 

21 

33 

68 

172 

30 

44 

93 

179 

32 

46 

116 

191 

41 

56 

143 

214 

41 

57 

168 

238 

50 

66 

198 

270 

50 

67 

226 

298 

59 

77 

- 

- 

59 

78 

- 

- 

1 

h 

32 

38 

40 

187 

33 

39 

75 

177 

43 

49 

99 

186 

43 

50 

122 

197 

53 

60 

149 

220 

53 

60 

173 

243 

62 

70 

203 

275 

63 

70 

- 

- 

72 

80 

- 

- 

72 

80 

_ 

— 

ic 

Inventory  Level.  H  Policy 


H  <  Hl 

Buy  -  H 

H>  H2 

Dispose  H  -  K, 

**  Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n, 
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TABLE  3-A  Cont'd 

Policies  Resulting  from  Model  Runs* 


Run  509 


Period** 

h 

\ 

10 

10 

32 

40 

187 

9 

14 

38 

72 

176 

8 

22 

49 

99 

184 

7 

22 

50 

122 

196 

6 

31 

70 

151 

222 

5 

31 

94 

178 

247 

4 

39 

105 

209 

280 

3 

40 

115 

- 

- 

2 

48 

126 

- 

- 

1 

48 

135 

- 

- 

Run  511 

Period** 

\ 

h 

h 

10 

19 

36 

- 

- 

9 

22 

40 

- 

- 

8 

30 

50 

- 

- 

7 

30 

63 

- 

- 

6 

40 

91 

- 

- 

5 

41 

104 

- 

- 

4 

49 

116 

- 

- 

3 

50 

126 

- 

- 

2 

58 

138 

- 

- 

1 

59 

148 

- 

- 

Run  510 


Period** 

h 

h 

h 

h 

10 

28 

39 

40 

187 

9 

29 

40 

76 

178 

8 

38 

50 

107 

191 

7 

38 

50 

129 

203 

6 

48 

70 

157 

229 

5 

48 

109 

185 

253 

4 

57 

120 

216 

287 

3 

58 

130 

- 

- 

2 

67 

140 

- 

- 

1 

67 

150 

- 

- 

Run  512 

Period** 

1 1 

h 

10 

29 

39 

40 

41 

9 

30 

40 

80 

81 

8 

39 

50 

122 

123 

7 

49 

63 

156 

157 

6 

48 

96 

192 

193 

5 

49 

111 

227 

228 

4 

58 

122 

268 

269 

3 

59 

132 

- 

- 

2 

67 

143 

- 

- 

1 

68 

153 

- 

- 

Inventory  Level.  H  Policy 


H  < 

Buy  Kx  -  H 

H  >  H2 

Dispose  H  -  K, 

**  Real  time  periods;  i.e.  period  n  +  1  is  later  than  period  n 


APPENDIX  B 

FORTRAN  PROGRAMS  FOR  INVENTORY  MODELS 
DEVELOPED  IN  CHAPTER  4  AND  6 


The  listings  of  the  FORTRAN  programs  used  for  calcula¬ 
ting  optimal  policies  in  Chapters  4  and  6  are  presented 
in  this  appendix.  The  program  decks  are  available  on 
request. 


B-2 


FORTRAN  Program  for  Dynamic  Procurement 
Disposal  Model  in  Chapter  IV. 
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Appendix  C 

Discussion  of  ASO  Procedures 

A.  Demand  ("Requirement")  Prediction 

Two  systems  of  demand  prediction  are  in  use  at  ASO:  The  Program 
Usage  Replenishment  System  (PURS)  and  the  Replenishment  Demand  Issue 
System  (RDIS).  Under  both  systems  total  demand  is  estimated  for  a  time 
period  equal  to  lead  time  (production  lead  plus  activity  lead  plus  ad¬ 
ministrative  lead)  plus  a  six  month  "safety  level"  time  plus  a  six  month 
"order  cycle".  Safety  level  and  order  cycle  do  not  vary  from  item  to 
item,  although  lead  time  does  somewhat.  The  first  two  numbers  are  es¬ 
tablished  by  fiat  and  are  independent  of  costs  associated  with  items 
(holding,  purchase,  stockout,  etc.);  so,  since  the  lead  time  cannot  be 
manipulated,  it  is  clear  that  there  is  no  room  for  optimization  under 
either  PURS  or  RDIS. 

Under  PURS,  past  usage  records  are  kept  in  terms  of  items  consumed 
by  aircraft  maintenance  cycle  (24C  hours  for  all  aircraft  except  jets, 
for  which  it  is  60  hours).  These  records  are  divided  into  a  maintenance 
rate  (items  consumed  per  flying  hour)  and  an  overhaul  requirement  (items 
consumed  in  overhaul  at  the  end  of  the  cycle).  These  records  are 
gathered  from  user  activities  in  the  field. 

From  CNO  and  BUWEPS  come  estimates  of  expected  flying  hours  and  over¬ 
haul  schedules  during  the  following  six  months  and  for  the  ninth  month 
from  payment.  A  technical  file  is  maintained  which  provides  information 
on  item  applications.  There  appears  to  be  a  very  slight  relationship 
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between  active  lead  time  and  the  figure  used  in  these  calculations. 

PURS  then  uses  the  usage  records,  expected  future  flying  hours,  and 
the  technical  file  to  produce  a  requirement  for  each  item. 

Two  major  weaknesses  are  inherent  in  PURS.  The  first  is  the  arbi¬ 
trariness  of  the  demand  period  considered  and  the  second  is  the  unre¬ 
liability  of  the  data.  Usage  records  are  poorly  kept  at  activities  and, 
further,  no  significant  correlations  have  been  found  either  between  pre¬ 
dicted  flying  hours  and  actual  flying  hours  or  between  actual  flying 
hours  and  demand. 

Under  the  second  demand  prediction  system,  RDIS,  the  assumption  is 
made  that  future  demand  will  be  the  same  as  past  demand.  This  system 
has  the  advantage  of  simplicity,  and  it  does  not  rely  on  as  much  un¬ 
reliable  data  as  does  PURS.  However,  its  demand  estimates  are  no  better 
than  those  generated  by  PURS,  and  it  has  less  scientific  appeal; 

B.  Requirement  Adjustment  for  Repairables 

For  repairable  items,  predicted  demand  is  reduced  by  predicted  re¬ 
turns  of  items  which  have  successfully  passed  the  maintenance  cycle. 

The  maintenance  cycle  is  defined  as  six  months  for  most  items,  though 
only  three  for  some  high  value  articles.  Time  spent  by  an  item  in  the 
cycle  is  referred  to  as  "turnaround  time".  The  cycle  consists  of  three 
states:  a)  the  item  is  removed  from  an  aircraft  and  screened  for  re- 
pairability.  If  it  looks  satisfactory,  then  b)  it  goes  to  the  main¬ 
tenance  depot  where  it  is  screened  again.  If  it  is  still  alright,  it 
is  entered  in  the  stock  records  and  c)  goes  to  the  overhaul  shop. 
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There  it  gets  the  final  stamp  of  approval,  is  repaired,  and  is  put  on 
the  shelf. 

The  "recovery  rate"  is  defined  as  the  fraction  of  items  removed  from 
aircraft  which  return  from  the  maintenance  cycle.  The  "RB  recovery  rate" 
is  the  fraction  of  those  which  are  not  discarded  at  removal  which  return 
from  stages  b)  and  c)  of  the  cycle. 

The  order  requirement  for  repairables  (0)  is  calculated  as  gross  ex¬ 
pected  demand  (D)  minus  expected  item  removals  (R^)  times  recovery  rate 
(r)  minus  damaged  equipment  (d)  reported  on  hand  times  RB  (RB)  recovery 
rate  minus  amount  (H)  on  hand  and  on  order. 

0«D-Rixr-dxRB-H 

The  time  period  considered  for  demand  is  as  defined  for  PURS  and 
RDIS,  the  sum  of  lead  time,  safety  level  and  order  cycle.  The  removal 
period  begins  a  turn-around  time  earlier  and  ends  a  turnaround  time  plus 
a  safety  level  earlier. 

C.  Life  of  Type  Items 

At  one  time  the  Navy  estimated  its  requirement  for  some  items  by 
"calculating"  total  demand  overall  time  and  reducing  this  by  estimated 
recoveries.  This  method  was  usually  used  for  high  recovery  rate  items. 
At  present,  however,  they  buy  only  eighteen  months  worth  of  such  items. 
This  generally  amounts  to  the  same  thing,  when  safety  factors  have  been 
taken  into  account. 

ASO  has  a  formula  which  it  uses  to  decide  whether  to  dispose  of  some 


life-of-type  aircraft  parts.  The  formula  produces  the  expected  number 
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of  maintenance  cycles  in  the  remaining  life  of  the  aircraft  to  which  the 
part  is  applied;  this  number  multiplied  by  the  appropriate  usage  rate 
produces  a  desired  number  of  parts  on  hand.  If  this  is  less  than  the 
actual  number  on  hand,  then  ASO  will  dispose.  The  formula  is: 


Number  of  maintenance  cycles  in  life  of  item 


where  N  is  number  of  aircraft,  h  is  average  monthly  flying  hours  per 
aircraft,  a  is  aircraft  survival  rate,  e  is  system  recovery  fraction, 
x  is  number  of  months  of  remaining  life  averaged  over  all  aircraft  cal¬ 
culated  from  service  tours  remaining,  k  is  number  of  months  turnaround 
time  plus  safety  level,  g  is  fraction  of  aircraft  operating,  and  j  is  a 
safety  factor  which  ASO  has  introduced  as  a  hedge  against  the  possible  un¬ 
reliability  of  this  formula. 


Economic  Order  Quantity 


A  procedure  for  calculating  an  "economic  order  quantity"  (EOQ)  is 
applied  to  some  consumable  items.  It  produces  an  operating  level  of 
an  ite.::  to  which  lead  time,  etc.,  are  then  added.  The  EOQ  is  required  to 
lie  between  the  operating  level  calculated  on  the  basis  of  estimated  demand 
over  the  time  period  mentioned  above  and  the  life  of  type  estimate  of  the 
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item.  It  is  designed  to  correct  egeinst  uneconomic elly  smell  orders, 
but  cannot  prevent  the  more  usual  uneconomical ly  large  orders.  The 
EOQ  procedure  also  provides  a  scientific  method  for  computing  a  safety 
level  as  follows: 

a)  One  computes  D^L  +  L  D  where  D  is  expected 

demand  and  L  is  expected  lead  time.  The  first  L 
and  the  second  D  in  the  equation  are  approximations 
to  standard  deviations  of  demand  and  lead  time. 

b)  One  computes  R  -  (Q  x  I  x  C)/(T  x  H  x  E)  where  Q 

is  the  economic  order  quantity,  I  is  a  discount  fact¬ 
or  which  includes  Interest,  obsolescence,  and  holding 
cost,  C  is  the  unit  price,  T  is  estimated  annual  de¬ 
mand,  H  is  the  shortage  cost,  and  E  is  "military 
essentiality." 

c)  One  uses  I-R  and  a  table  of  normal  variates  to  find 
K,  the  number  of  multiples  of  standard  deviation. 

d)  One  declares  that  the  safety  level  is  K  times  o^L, 

This  is  assumed  to  be  in  units  of  quantity,  dimensional 
analysis  to  the  contrary  notwithstanding. 

E.  Critique  of  These  Procedures 

A  major  handicap  in  the  Navy's  thinking  about  supply  procedures 
lias  in  the  depiction  of  expected  demand  estimates,  based  on  suspect 
data  and  an  arbitrary  time  frame,  as  "requirements."  Once  a  course  of 
action  having  the  repercussions  of  this  one  has  been  called  a  requirement 
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to  be  met  at  all  cost,  then  no  quantity  of  peripheral  models  and  tech¬ 
niques  will  introduce  a  significant  increase  in  supply  efficiency.  If 
the  Navy  has  any  desire  to  reduce  the  huge  cost  of  its  supply  operation 
it  must  have  a  system  which  considers  simultaneously  demand,  purchase 
costs,  holding  costs,  and  outage  costs,  and  which  neither  equates  quan¬ 
tities  with  times  nor  substitutes  mean  demands  for  distributions  of 
demands . 

If  the  concept  of  "safety  'level"  had  any  operational  meaning,  then 
it  clearly  should  vary  from  item  to  item.  In  fact,  however,  it  does  not 
have  meaning.  For  a  given  state  of  the  inventory  syBtem,  there  exists 
a  best  quantity  of  each  item  to  order.  Other  things  remaining  constant, 
this  best  quantity  will  vary  with  the  production  set-up  cost  of  the 
item.  If  there  were  Such  a  thing  as  a  safety  level,  then,  it  would  have 
to  be  something  which  is  a  function  of  set-up  cost.  This  is  clearly  a 
contradiction  in  terms. 

F.  The  Data  System 

There  are  two  major  weaknesses  in  the  ASO  data  collection  and  main¬ 
tenance  procedures.  The  first  is  the  unreliability  of  what  is  collected 
and  the  second  is  the  choice  of  what  ought  to  be  collected.  The  un¬ 
reliable  data  problem  is  one  which  ASO  can  hardly  solve  by  itself,  but 
which  must  be  faced  by  the  Navy  as  a  whole;  however,  something  can  be 
done  about  deciding  what  data  would  be  useful  to  have. 

Certainly  summaries  of  the  obsolescence  experience  of  items  should 
be  kept--e.g.,  for  each  category  of  item,  the  average  lifespan  of  its 
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members.  Similarly,  summaries  of  demand  history  should  be  maintained. 
Stockout  reports  should  be  input  to  the  data  collection  system.  This 
sort  of  information  will  be  required  by  any  effective  Inventory  control 
model,  and,  together  with  such  a  model,  could  enable  a  truly  efficient 
operation  to  be  carried  on. 
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